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ABSTRACT 


This study investigates the influence of a curvilinear fiber format on load car- 
rying capacity of a layered fiber-reinforced plate with a centrally located hole. 
A curvilinear fiber format is descriptive of layers in a laminate having fibers 
which are aligned with the principal stress directions in those layers. Lami- 
nates of five curvilinear fiber format designs and four straightline fiber format 
designs are considered. A quasi-isotropic laminate having a straightline fiber 
format is used to define a baseline design for comparison with the other lami- 
nate designs. Four different plate geometries are considered and differen- 
tiated by two values of hole diameter/plate width equal to 1/6 and 1/3, and two 
values of plate length/plate width equal to 2 and 1. With the plates under 
uniaxial tensile loading on two opposing edges, alignment of fibers in the 
curvilinear layers with the principal stress directions is determined analytically 
by an iteration procedure. In-plane tensile load capacity is computed for all 
of the laminate designs using a finite element analysis method. A maximum 
strain failure criterion and the Tsai-Wu failure criterion are applied to deter- 
mine failure loads and failure modes. Resistance to buckling of the laminate 
designs to uniaxial compressive loading is analyzed using the commercial 
code Engineering Analysis Language. Results indicate that the curvilinear fi- 
ber format laminates have higher in-plane tensile load capacity and compara- 
ble buckling resistance relative to the straightline fiber format laminates. 
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CHAPTER I 


INTRODUCTION 

Advanced fiber-reinforced composite materials were in- 
troduced roughly two decades ago as a high performance 
structural material with enormous possibilities for increas- 
ing structural efficiency. Initial research on composite 
materials centered on producing the material at a reasonable 
cost and fabricating simple structural shapes from the ma- 
terials. Additional research was also conducted with regard 
to analyzing and predicting how those simple structural 
shapes would respond to loads . 

The next area of research focused on the complex and 
sometimes unpredictable failure mechanisms of composite ma- 
terials. Many variables affect the structural response of 
composite materials, and since the behavior of metals had 
been well studied, many fiber-reinforced components were de- 
signed to behave somewhat like metals. This approach was 
acceptable because fabrication of the components was of 
major concern. The metal-like construction of composites 
was accepted because the lower density of composite mater- 
ials produced a weight savings in the components. However, 
use of composite materials in this manner did not take ad- 
vantage of the major attributes of composite materials which 
allow for tailoring the composite material to fit the appli- 
cation. The strength and stiffness of a composite material 
can be aligned in directions which correspond to the loads 
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applied to a structural component. 

Military aircraft and both military and civilian space- 
craft have been designed and used successfully with advanced 
fiber-reinforced composite materials as part of the overall 
design. Fiber-reinforced composite materials are now being 
studied for use in business aircraft, and to some extent, 
commercial aircraft. Evern though composite materials have 
seen increased usage , the potential advantages of advanced 
fiber-reinforced composite materials have not been fully 
realized, or even investigated. The development and use of 
composite materials in structures has followed a conserva- 
tive approach. This conservative approach has resulted in 
some of the detrimental characteristics of composites erod- 
ing some of its advantages. As a result, in many cases only 
a marginal gain in efficiency is realized when a composite 
material is used instead of a metal. For example, a quasi- 
isotropic laminate has stiffness characteristics which are 
similar to aluminum and has similar in-plane load capacity. 
However, unlike aluminum, a quasi-isotropic laminate may 
delaminate, while aluminum does not have that problem. A 
quasi-isotropic laminate is susceptible to environmental 
degradation, while aluminum has less of a problem with this. 
The manufacture and repair of composite materials is less 
fully understood than the manufacture and repair of alumi- 
num. The addition of more weight to a composite structure 
might alleviate these problems, but then the advantages of 
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using composite materials begin to decrease. Gains made in 
structural efficiency have been moderate while it seems the 
gains could be far greater. To make gains in structural 
efficiency by a factor of 2 or 3, the overly conservative 
utilization of a composite material must be disgarded and 
more innovative approaches must be developed. This parti- 
cular study examines one departure from the conservative use 
of composite materials and develops an approach which may 
lead to step increases in structural performance. 

Conventional design philosophies for fiber-reinforced 
composite structures are based on the idea of using multiple 
layers of fibers embedded in a matrix, the fibers in each 
layer being straight and parallel to each other, and aligned 
in a particular direction. Though each layer may have its 
own unique fiber orientation, the idea of allowing the fiber 
orientation within a layer to vary from point to point has 
not been considered seriously. The specific issue in this 
study is to use fiber reinforcement in such a way that the 
direction of the fibers, or at least some of the fibers, is 
a function of spatial position in the structure. This will 
be referred to herein as a curvilinear fiber format. The 
curvilinear format has the following advantage: A structur- 
al component may contain a geometric discontinuity, such as 
a hole, which interrupts fiber continuity in all of the lay- 
ers and causes a concentration and realignment of stress. 
If the fibers are straight and parallel to each other, a 
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large number of fibers must be broken at this geometric dis- 
continuity . The efficient way to use the fibers would be 
not to break the fiber continuity and thus use the fibers to 
greater advantage near the geometric discontinuty . The 
fibers should "flow" continuously around the discontinuity 
and be oriented in such a manner as to transmit the load 
efficiently around the (discontinuity. Many issues must be 
studied to make this idea viable, the most important issue 
being fabrication. However, the availability of raw fiber, 
the increased power and flexibility of robotics, and new 
matrix development do provide promise for fabrication of 
components on a fiber-by-fiber basis with something other 
than a straightline format.^ 

In considering the issue involved, this study focuses 
on the application of the concept to a particular problem. 
The particular problem studied here is a prime candidate for 
deviating from the straightline fiber format. The particu- 
lar problem is a thin plate which contains a centrally lo- 
cated circular hole and is loaded uniaxially in its plane. 
The problem has been studied hundreds of times in the con- 
text of using isotropic and composite materials. However, 
it is an important problem. Aircraft structures contain 
many holes for access and fabrication, and commercial air 
transports contain numerous windows. The central question 
is: If the fibers are not to be used in a straightline 
fiber format, how should the orientation of the fibers vary 
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from point to point? The basic philosophy used here to de- 
sign the plate will be to align the fibers with the princi- 
pal stress directions. More specifically, in those layers 
which are allowed to have variable fiber orientation, the 
fibers will be aligned with the principal stress directions 
in those layers. The method focuses on the use of an itera- 
tion scheme to align the fibers with the principal stress 
directions. Several different stacking sequences of layers 
and four different geometries are considered to see how the 
curvilinear fiber format will compare with the straightline 
fiber format. The maximum strain and Tsai-Wu failure cri- 
teria are used to determine the comparative strengths and 
weaknesses of the various designs. Finally, the buckling 
strengths of the curvilinear designs are computed to deter- 
mine if the curvilinear designs alter buckling strength. 

The next chapter begins to address these issues by pre- 
senting a brief overview of classical laminated plate 
theory, casting these results within the framework of varia- 
tional methods, and deriving the governing nonlinear partial 
differential equations and associated boundary conditions. 
Chapter III presents the specific plate geometries to be 
studied, and presents the governing equations developed in 
Chapter II in the context of those plates subjected to in- 
plane uniaxial tensile loading. The equations, by the 
nature of the problem, are linear partial differential 
equations with coefficients that are a function of the spa- 
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tial coordinates. The material properties vary from point 
to point and these material properties modify the differen- 
tial equations. Because of the variable coefficients, an 
approximate numerical scheme is necessary to obtain numbers. 
Here the finite element method is used, and details of the 
particular finite element used are presented in Chapter III. 
Since the finite element' method is well known, the discus- 
sion of the finite element approach is brief. The boundary 
conditions used in the tensile load analysis and the use of 
a one-quarter plate analysis are also discussed. Most 
importantly, the iteration scheme used to align the fiber 
directions with the principal stress directions is discuss- 
ed. Finally, Chapter III includes an overview and discus- 
sion of the failure thoorles used. 

Chapter IV focuses on the buckling analysis. The 
method of adjacent equilibrium is used to derive the equa- 
tions which govern buckling. The buckling equations are 
based on the equations derived in Chapter II. Again these 
equations are partial differential equations with variable 
coefficients and solutions are obtained using the finite 
element method. Since this particular finite element analy- 
sis is based on use of the commercial code Engineering Anal- 
ysis Language (EAL), the discussion centers on the charac- 
teristics of the element used. The boundary conditions used 
in the buckling analysis are presented along with reasons 
for analyzing the entire plate instead of only one-quarter 
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of the plate. 

Chapter V presents the results for the tensile load 
problem. The failure loads and failure modes for 36 diff- 
erent plate designs are presented. Since the point of the 
work is improved performance with the curvilinear designs, 
gains in using this design are discussed in some detail. 
Through the use of contour plots, the basis for the improve- 
ment in terms of the redistribution of stress is discussed. 
Convergence and accuracy of the stress analysis for the lam- 
inates with the curvilinear layers is also addressed in 
Chapter V . 

Chapter VI centers on the buckling resistance of the 36 
different plate designs. As mentioned before, the purpose 
of the buckling analysis was not to design for improved 
buckling. Rather the buckling analysis was conducted to de- 
termine the influence of the curvilinear fiber format on the 
buckling capacity. 

Finally, the last chapter. Chapter VII, includes a sum- 
mary of the findings, discusses the limitations of the re- 
sults, and makes recommendations for future research. 
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CHAPTER II 


DERIVATION OF EQUATIONS GOVERNING 
GENERALLY ORTHOTROPIC 
THIN LAMINATED PLATES 

General Considerations 

Consider a thin laminated fiber-reinforced plate as 
shown in Figure 2.1. The x-y coordinate plane corresponds 
to the middle surface of the plate. The z-axis is perpen- 
dicular to the middle surface and z is measured relative to 
the middle surface. The total thickness of the plate is 
denoted by H . 

The plate studied here is governed by the following set 
of assumptions based on classical thin-plate theory. 

(1) The thickness of the plate is uniform and small compar- 
ed to the other dimensions of the plate. 

(2) The Kirchhoff hypothesis is assumed to be valid. The 
hypothesis states that at any point on the middle surface, 
line elements normal to the undeformed middle surface remain 
straight, normal to the surface, and do not extend or con- 
tract during the deformation of the plate. As a result, the 
transverse shear strain components are exactly zero. This 
is illustrated in Figure 2.2 where the orthogonality of the 
line element relative to the middle surface is shown. 

(3) The strains are small compared to unity and the rota- 
tions about the x-axis and the y-axis are moderately small. 
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The shrains aire sroalleir hhan the rotations . The rotations 
about the z-axis are negligibly small. Essentially it is 
assumed that the von Karman thin-plate theory is valid for 
this problem, specifically the buckling problem. 



Figure 2.1 Thin Plate 


9 





Figure 2.2 Kirchhoff Assumption 
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(4) The plate is constructed of multiple layers of fiber- 
reinforced material, each layer being bonded perfectly to 
tho adjacent layer. Each layer is of thickness h and has 
its own fiber orientation. 

(5) The plate is in a state of plane stress, i.e., at every 

point in the plate, a, = t,, = = 0. 

(6) The stiffness properties are uniform throughout the 
thickness of each layer, but the properties are a function 
of locations x and y in each layer. Each layer can also 
have different stiffness properties. 

(7) Each layer exhibits linear elastic behavior. 

(8) The material in each layer is orthotropic. 

(9) There are no body forces acting in the plate. 

Detailed Considerations For Layered Fiber-Reinforced Plates 
In terms of plate displacements, the Kirchhoff assump- 
tion can be written as: 


u(x,y, 2 ) 

= u0(x,y) 

z (wO(x,y) ) 

1 * 

(2.1a) 

v(x,y,z) 

= v0(x,y) 

z (wO(x.y)) 

I J 

(2.1b) 

w(x,y,z) 

= w0(x,y) 


(2.1c) 


where u, v, and w denote the displacements at any point in 
the plate, and u° , vO , and w® denote the displacements of 
the middle surface of the plate. The quantity wO, is the 

t " 

rotation of the normal to the middle surface about the y- 

axis, as illustrated in Figure 2.2, and the quantity wO is 

» y 

the rotation of the normal to the middle surface about the 
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x-axis , also illustrated in Figure 2.2. 

The strain-displacement relations, including the 
effects of moderate rotations about the x and y axes are: 


(x,y,z) 

= (u(x,y,z)) , 

+ 

J^CCwCx.y.z)) „]2 

(2.2a) 

€y (x,y,z) 

= (v(x.y,z)) 

» 7 

+ 

»$[(w(x.y,z)) T2 

> 7 

(2.2b) 

(x,y, z) 

= ( u ( X , y , z ) ) , 

1 7 


(v(x,y , z) ) ^ 



+ C (m(x, 7 ,z)) (w(x,y,z)) ] (2.2c) 

I * i 7 

where are the strains at any point in the plate. 

Substituting the Kirchhoff displacements into the strain- 
displacement relations results in: 


€ = uo 

* » X 


+ 

H [wO 

(2.3a) 

€ = vO 

7 . y 

- “ "°.T7 

+ 

^ [mO „]2 

1 7 

(2.3b) 

r = 

*y . y 

+ yO 

1 X 

2 

Z W® ,v 

> • 7 

.x«°.y (2.3c) 

Written in more compact notation: 



6^ (x,y ,z) = 

€0, (x.y) 

+ " 

z (x.y) 

(2.4a) 

Gy (x.y.z) = 

€0y (x.y) 


z K®y (x.y) 

(2.4b) 

r^y(x,y,z) = 

r®^y (x,y) 

+ 

2 KOjjy(x.y) 

(2.4c) 

where the strains 

6° of the 

middle surface 

of the plate 

are : 





€.0 = uO 

X .X 

+ [wO 

f * 

]2 


(2.5a) 

€0 z yO 

7 .7 

+ CwO _ 

1 7 

]2 


(2.5b) 

o 

II 

c 

o 

+ yO + 

1 X 

W 

0 uO 

1 X , y 

(2.5c) 

and the curvatures 

K° of the 

middle surface 

of the plate 

are : 





= - wO 

X , , 

X X 



(2.6a) 
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K°y = - W®,yy (2.6b) 

K°xy = - 2 (2.6c) 

A single layer within the plate is often referred to as 
a lamina. A lamina consists of many thousands of strong and 
stiff fibers embedded in a softer and weaker material . This 
second material is referred to as a matrix. The fibers 
carry most of the load while the matrix material transfers 
the load between the fibers and holds the fibers together. 
A laminate is a group of laminae or layers that, as stated 
before, are bonded perfectly together. The planes of the 

laminae are parallel to each other. This is shown in Figure 

2.3. In Figure 2.3 the fiber directions of the various 
laminae are illustrated, and also as stated before, each 
lamina can have its fibers oriented in an arbitrary 

direction . 

Because each layer can have its own fiber orientation, 
it is convenient to use a global or laminate coordinate 
system for the analysis. The global system is referred to 
here as the x-y-z coordinate system while the lamina system 
is denoted as the 1-2-3 coordinate system. The relationship 
between the 1-2-3 lamina coordinate system and the x-y-z 
laminate coordinate system is shown in Figure 2.4. The 
1-axis corresponds to the fiber direction. The angle 0 
defines the orientation of the 1-axis with respect to the 
x-axis . 

One of the most important aspects of laminated fiber- 
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+ 9 : Positive Rotation Of Coordinate Axes 


Figure 2.4 

Relationship Between 1-2-3 Lamina Coordinate 
System and x-y-z Laminate Coordinate System 


reinforced composite plates, and one which distinguishes 
them from traditional isotropic plates, is the stress-strain 


behavior. As mentioned on page 11, it is assumed that the 
material in each layer is orthotropic and exhibits linear 
elastic behavior. Referring to Figure 2.5, which shows the 
stresses acting on an elemental unit cube in the 1-2-3 
coordinate system, the s'tress-strain behavior of an ortho- 
tropic material in the 1-2-3 coordinate system is given by: 


23 


3 1 


1 2 


’ll 


I 2 


I 3 


0 

0 

0 
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'2 2 


'2 3 


0 

0 

0 


'1 3 


'23 


'3 3 


0 

0 

0 


4 4 


5 5 


6 6 



■ ^l 


^2 


^3 


^2 3 


Tai 


^1 2 


(2.7) 


where the C 's are the stiffness coefficients of the mater- 
ial. 


The stress-strain relations of Eqn . 2.7 reduce to the 
following form for an orthotropic lamina in a state of plane 
stress . 


b 


Qi 1 2 0 


1 

m 

<^2 


^12 ®22 ^ 


^2 

1 


_ 0 0 Qee 


r ,2 


( 2 . 8 ) 


The Q's are 
coordinate system. 


the reduced stiffnesses in the 1-2 
The Q's can be defined in terms of the 
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1 


Figure 2.5 

Stresses Acting On An Elemental Unit Cube 
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engineering constants , E2 

E, 


u 


1 2 




2 1 


, and G^2 • 


Qii = 


1 2 


22 


1 

- ^^ 12^^21 



^ 2^12 


1 

^ 12^21 

1 - ^^^ 2^^21 


E2 


1 

- ^ 12^21 ' 



( 2 . 9 a) 


( 2 . 9 b) 


( 2 . 9 c) 


86 " *^12 


( 2 . 9 d) 

To completely describe the response of an orthotropic 
material in a state of plane stress, four independent mater- 
ial properties are required. These material properties are 
the elastic modulus in the fiber direction , the elastic 
modulus perpendicular to the fibers and in the plane of the 


of the lamina 


, the in-plane shear modulus 


1 2 


, and 


Poisson's ratio Poisson's ratio 421 

terrained from the following reciprocal relation: 


4 


1 2 


4 


2 1 


(2.10) 


In general, Poisson's ratio 4. . measures the contraction 

* w 

in the j-direction due to a stress applied in the i-direc- 
tion. 

The reduced stiffnesses in the x-y coordinate system 
are determined by transforming the reduced stiffnesses in 
the 1-2 coordinate system into the x-y coordinate system. 
The transformation matrix [ T ] needed for this is: 
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cos 2 0 


sin2 0 


2 sin 0 cos 9 


sin2 0 cos2 0 -2 sin 9 cos 0 

-sin 0 cos 0 sin 0 cos 9 cos2 0 - sin2 0 


( 2 . 11 ) 


where 0 denotes the angle from the x-axis to the 1-axis 
for a positive rotation in a Cartesian coordinate system. 
Given values of stress or strain with respect to the x-y 
coordinate system, these may be transformed to the 1-2 
coordinate system according to the following: 


1 

Q 

1 


■ 


= [ T ] 


Ti2 


H 

M 


(W 

1 


■ 

CM 

= [ T ] 


^12 



2 


2 


Conversely, given stress or strain in the 1-2 coordinate 
system, these may be transformed to the x-y coordinate 
system according to the following: 


1 

Q 

M 

‘ 




1 

m 

M 


1 


= [ T ]-i 

^2 



= [ T ]-i 

^2 



^12 


y 


ri2 




2 


2 


The reduced stiffnesses in the x-y coordinate system 
for a lamina, given the reduced stiffnesses in the 1-2 
coordinate system, can be determined using the following 
equation: (Jones [1], p. 50). 
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[ Q ] 


[ T ]-i [ Q ] C T ]-T 


( 2 . 12 ) 


where : 

[ Q 3 = reduced stiffnesses in the x-y coordinate system 

[ T 3"i = inverse of transformation matrix [ T ] 

[ Q 3 = reduced stiffnesses in the 1-2 coordinate system 

[ T 3"'^ = transpose of inverse of matrix [ T 3 

The stress-strain relations in the x-y coordinate system for 
an orthotropic lamina in a state of plane stress then 
become : 




(2.13) 


It is important to note that in the x-y coordinate 
system there is coupling between the shear strain and the 
normal stresses, and coupling between the shear stress and 
the normal strains. The lamina appears to be anisotropic in 
x-y coordinates. However, the lamina is orthotropic because 
there are still only four independent material properties. 
The lamina is designated as generally orthotropic and can be 
characterized by equation 2.13. (Jones [13. P. 51). 
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Force And Moment Resultants 


The z-location of a layer in a laminate is shown in 

< 

Figure 2.6. In particular, the k-th lamina is bounded by 
the coordinates and When a plate is subjected to 

loads, stresses develop in each layer. To aid in deriving 
the governing differential equations, these stresses are 
lumped into equivalent forces and moments, commonly referred 
to as resultants . The force and moment resultants acting on 
an element of the plate are shown in Figure 2.7. These 
resultants have units of force per unit length or moment per 
unit length in the x and y directions. and are in- 
plane normal force resultants, is the in-plane shear 

force resultant, M^^ and M^ are bending moment resultants, 
and is the twisting moment resultant. 

The force and moment resultants are related to the 
stresses by the following equations: 


= 1 ^x (2.14a) 

. -hh 

Ny = I Oy dZ (2.14b) 

-HH 

N., = I -fxydz (2.14c) 

->iH 


= I a, z dz (2 . 14d) 

-KH 

My = I Oy z dz (2.14e) 

-HH 

4-y,H 

M^y = J T^yZ dz (2. 14f ) 

-HH 
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Figure 2. 


6 Geometry For Layers In A Laminate 
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Substituting equation 2.13 into equations 2.14a-f 

results 

in: 


+ HH 


+ HH 

+>$H 



Nx 

= 

J Qii^x dz 

+ 

J Q^ 2 €, dz + J 

f ^16 r,, dz 


(2.15a) 



-HH 


-HH 













= 

J 

+ 

I ^ 2 ^y dz + 

[ ^26 r,y dz 


(2.15b) 



-SH 


-hh 

->5H 





+ HH 


+ HH 





= 

1 9>.«. dz 

+ 

I 926^, dz * J 

Q«e r dz 

^8 6 X y 


(2.15c) 



-SH 


-HH 

-^H 





+ HH 


♦ HH 




Mx 

- 

J ^ 11 ^x 2 dz 

+ 

J Qi 2 €yZ dz + 

J ®18 ^xy^ 

dz 

(2 . 15d) 



‘ -HH 


-HH 

-}$H 





+ HH 


♦ HH 

+^H 




= 

J dz 

+ 

I ® 22 ^y 2 dz + 

J ^26 

dz 

(2 . 15e) 



-HH 


-HH 

-^H 








+^H 




z 

J Qie^x^ dz 

+ 

1 Q28^y2 dz + 

J ®66 ^xy^ 

dz 

(2 . 15f ) 



-HH 


-HH 





The strains ^xy expressed in 
terms of the middle surface strains ^®x * ^°y » ^°*y 
and middle surface curvatures K<>, , KO^ , K® using 
equations 2.4a-c. These expressions can be substituted into 
equations 2.15a-f, Since the middle surface laminate quan- 
tities are not functions of z , equations 2.15a-f become; 
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nr fS 

OF POOR QUAUTV 


N. 


+HH +KH 

6”. J + K“. 1 2 02 


-HH 


HH 


+ HH 


+ HH 


^°y J '^12 ^ K®y J ®12 ^ dz 


-HH 


■HH 


♦HH :. 

* ro,^ I dz + K0„ J z dz 

-HH 


+ KH 

Q 

-HH 


( 2 . 16 a) 


+HH 


+ HH 


N. 


J ®12 ■*■ 1 ^12 ’ 


-HH 


-HH 


-►HH 




N 


xy 



1 0=2 

dz * KO, J 

®22 2 

dz 


-HH 


HH 








1 

Qjg dz + 

KO f 

*y J 

®26 



HH 


HH 

+ HH 

+54H 



€0, 1 Q 

16 

+ KO, j z dz 


-HH 


-HH 








+ €0^ 

J «ae 

dz + KOy 

j Oas 

z dz 


-HH 


-HH 




♦ HH 


♦ HH 

+ 

ro 

xy 

J Qgg dz + 

KO 

*y 

1 ®6 



-HH 


-HH 


z dz ( 2 . 16 b) 
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z dz + 


M. 


+ HH 

1 

-HH 


+ HH 

KO, J Q.. dz 

-HH 


'I’HH 




+ eOy J Qj 2 Z dz + KOy J 22 dz 


-HH 


-54H 


♦ HH 




+ ro__ J Qjg z dz + Ko^y J Qjg z2 dz (2.16d) 


»y 


-HH 


-54H 


M. 


+HH +HH 

J ®X 2 z dz + KO, J Qj 2 z 2 dz 

-HH -»4H 


M 


* y 



♦ HH 





> 60^ 

1^2 2 

z dz 

^ K% 1 

^22 

Z2 dz 


-HH 


- 

HH 




+hh 



♦ HH 


r^y 1 

^2 S 

z dz + 

KO 

X y 

1 ®38 



HH 



-HH 

+ HH 


♦ HH 



€0, J Q 

16 2 dz + 

J Qi6 

dz 

-HH 



-HH 









+ €0^ 

1 ^26 

z dz 

+ KOy 

1 ^26 

z2 dz 


-HH 



-HH 




>HH 



♦ HH 

+ 

ro 

»y 

J «Z8 

z dz + 

KO 

*y 

J «66 



-HH 



-HH 


z 2 dz ( 2 . 1 6e ) 
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Each layer can have its own unique fiber orientation, 
so the Q's may be different for each layer. However, the 
Q's are constant within a single layer, so the integrals in 
equations 2.16a-f can be converted to sums over the layers. 
The sums are defined to be: 


‘1 1 


12 


♦ HH 


J dz 


-HH 
+ >4H 


I Qi2 ^2 


-HH 


^ 1 ^^k ■ ^k- I ^ 

k = 1 


^ 2 ^ k ^^k ~ ^k- I ^ 

k= 1 


^16 


'2 2 




J Qi 6 dZ 


-HH 
+ >4H 


J ®22 dz 


^ 6 ^ k ^ ^k - I ^ 

k= I 


^ ^^2 2 ^ k ^ ^k ^k - 1 ^ 

k = 1 


'28 




I ^26 dz 


-SH 


N _ 

^ ^^26 ^ k ^ ^k “ ^k * 1 ^ 

k = 1 


'66 




I ^86 


-SH 


5 \ ^^k ^k- 1 ^ 

k= 1 


B 


1 1 




J* Q ^ J z d z 


-SH 


N - (Z^W ■ Z2 ) 

I (Qii)k ; 

k= 1 ? 


B 


1 2 


+ HH 


J Qi 2 2 dZ 


-HH 


H - (Z2 - z2 ) 

2 (Qi2)k ; 

k= 1 9 


(2.17a) 


(2.17b) 


(2.17c) 


(2. 17d) 


(2.17e) 


(2. 17f ) 


(2 . 18a) 


(2.18b) 
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<22^ - 


M _ (Z2 - 22 ) 

2 (Q22)^ ^ 

K=1 2 


N ■ _ (Z2 - Z2 ) 

Z (Q2«), ^ 

»t=l 2 


^-1 ^^88 


(22^ - Z2^,p 


J (Qll)k — 

k= I 3 


- 22k-1> 


u /^^2>k 
k= I 


- 23,.,) 


N _ (Z3, - Z3,_ , ) 

2 (Qi6>k 

k = I 3 


2 (Q,2)k 

k=l “ 


(23, - Z3^_^) 


^ (^26^k 

k= I 


(23, - Z\.i) 


^ (^6 e ^ k 


(23^ - Z3,_,) 
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The A's are the extentional stiffnesses for the laminate, 
the B's are the bending/extentional stiffnesses for the 
laminate, and the D's are the bending stiffnesses for the 
laminate. These stiffnesses, using index notation, are: 


A 


i J 


B 


i j 


D 


i J 


N _ 


2 (Qij)k 
k = 1 


2 


k=l 




3 


i.J = 1 , 2,6 

i,j = 1 . 2,6 

i.j = 1 . 2,6 


( 2 . 20 ) 

( 2 . 21 ) 

( 2 . 22 ) 


Equations 2.16a-f in terms of the A's, B's, and D's, now 
become 




All 


+ 

Bii 

KO 

X 

+ 

Al 2 

€0y + B,3 KO^ 










A 1 6 

KO,^ 

(2 . 23a) 


= 

Al 2 


+ 

Bi2 

KO 

X 

+ 

A2 2 

eo + B,, KO 
y 2 2 y 









+ 

A26 

ro + B KO 

^ xy °26 xy 

(2.23b) 


= 

Al 6 

e"x 


Bi 6 

KO 

X 

+ 

A26 

€0 + KO 

y 2 6 y 










Agg 

+ B,, KO,^ 

(2.23c) 

Mk 


Bii 

€®x 

+ 

Dll 

KO, 

+ 

Bi 2 

€0^ + D,3 KO^ 









+ 

Bie 

+ D,, K0,y 

(2.23d) 

«T 

= 

Bi 2 

^°X 


Di 2 

KO, 

+ 

®2 2 

^ D33 Ko^ 









+ 

^26 

ro + D,e KO 

xy 26 xy 

(2 . 23e) 


r 

Bi 6 



Die 


+ 

^26 

€0 + Do. KO 

y 2 6 y 









+ 

^66 

ro + D.. KO 

xy 66 xy 

(2.23f) 
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Using matrix notation, equations 2.23a-f are: 

^11 ^12 ^16 ®ll ®12 ®16 

^12 ^22 ^26 ®12 ®22 ®26 

^16 ^28 ^66 ®16 ® 2 S ®66 

(2.24) 

®ll ®12 ®18 ^11 ^12 ®18 

^12 ^22 1>22 ^ 2 B 

®18 ®26 ^66 ®16 ®26 ^66 ^° X 7 

J k 

The coupling stiffnesses are zero when the lami- 

nate is symmetrical . These terns imply that there is coup- 
ling between bending and extension. Subjecting a laminate 
which has coupling stiffnesses to extenslonal loads will 
cause the laminate to bend and/or twist in addition to 
stretching. Subjecting a laminate which has coupling stiff- 
nesses to moments will cause' the middle surface of the lam- 
inate to extend and shear, in addition to bending and/or 
twisting . 

The coupling stiffnesses are equal to zero when the 
laminate is symmetric in both geometry and material proper- 
ties about the middle surface of the laminate. (See Jones 
[1]. page 160) 

If the laminate is symmetric, equation 2.24 decouples 
and results in: 
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(2.26) 


1 

K 

x: 


^11 ^12 ^16 


1 

o 



^12 ^22 ®26 


KO^ 

Mxt 


^16 ^26 ®6G 




The presence of the and Ajg stiffnesses in equa- 

tion 2 . 25 shows that there is coupling between normal force 
resultants and shearing strain, and coupling between the 
shearing force resultant and normal strains. The presence 
of the D^g and Djg stiffnesses in equation 2.26 shows 
that there is coupling between bending moment resultants and 
twisting curvature, and coupling between the twisting moment 
resultant and normal curvatures. 

Nonlinear Equilibrium Equations And Associated Boundary 
Conditions 

Consider a portion of a laminated fiber-reinforced 
plate as shown in Figure 2.8. The x-y coordinate plane 
corresponds to the middle surface of the plate. The z-axis 
is perpendicular to the middle surface and, as stated be- 
fore, z is measured relative to the middle surface. The 
total thickness of the plate is denoted by H. The interior 
region of the plate is denoted by Q The boundary of the 
hole in the plate as seen from the z-direction is denoted by 
Fj . For our case, the boundary F^ does not have prescribed 
forces acting on it and no displacements are specified. The 
boundary around the plate as seen from the z-direction is 
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denoted by Fj . For simplicity F 2 is considered at this 
point to be a smooth curve. In reality, Fg is the boundary 
of a rectangular plate. The in-plane force resultants 
and are prescribed on Fj . acts in the direction 

normal to the boundary Fg and acts in the direction 

tangential to the boundary Fg • The normal load distribu- 
tion q(x,y) acts downward on the region Q perpendicular 
to the x-y plane of the plate. 

The total potential energy of the loaded plate x , is 
defined as the sum of the total strain energy U and the 
potential of the applied force V . 

X = U + V (2.27) 

The equilibrium equations and boundary conditions are deter- 
mined by setting the first variation of the total potential 
equal to zero. (Shames and Dym [2], pp. 401-408). That is 
6^^^x = 6^^\u + V) = 0 (2.28) 

The total strain energy is 



(2.29) 


(2.30) 


This integral is a point function of the upper limit . if 
this integral is independent of the path taken between 
the limits 0 and 6. . , meaning that t. . dG . , must be a 

perfect differential. When these conditions are met, then 
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du. = T^. de^. (2.31) 
Recall from page 11 the assumption that each layer exhibits 
linear elastic behavior, meaning that each stress component 
is linearly related to all of the strains by Hooke's Law 
which can be stated as: 



“ ^ijki 

^kl 


(2.32) 

where 

^ijkl " 

Cjtki 

Since T,j = Tj, 


and 

^ijkl “ 

J Ik 

since 



Substituting equation 2.32 into equation 2.29 results in 


“ = 111 [ J de.j ] dl (2.33) 

where the quantity inside the brackets is a perfect differ- 
ential and the coefficients C. are constants. Equation 
^.33 may now be expanded to take the form 

“ = IJJ [ '5 =tjkl «ki Ckiu ] di!. 

V 

jkl “ ^kl i j 

Thus, the following is established: 

“ = HI dl (2.34) 

V 

Substituting t. ^ for C. in equation 2.34 results 

in 


U = 



1 J 


1 J 


dY. 


(2.35) 


As stated on page 11, it is assumed that the plate is in a 
state of plane stress, so expanding equatipn 2.35 for this 
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case leads to 


0 = III [ * 2 (2'36) 

V 

Returning to contracted notation, 



and noting that = 14 r„ 

Substituting these relations into equation 2.36 results 
in the following equation for the total strain energy for 
linear elastic behavior 

“ = '4 III [ <^. r„ ] dl (2.37) 

V 


Using equations 2.8 and 2.9, the stress-strain rela- 
tions in the x-y coordinate system for linear elastic behav- 
ior, an orthotropic material, and a state of plane stress, 


are : 


a 


X 


a 


7 


1 "" u u 

^x X 




1 - U JJL 

y^y X 




(2.38a) 


(2.38b) 


"*”x y ^x y ^x y 

From equation 2,9b 


E u 

y*^xy 


1 UL U 

^xy^yx 


E jJL 

x^y X 

1 - u u 
^x y '^y : 


(2.38c) 


(2.39) 
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Substitute equations 2 . 38a-c into equation 2.37 to 


obtain 


U = 


III [ 

w ^ 1 - UL U 


E 


1 " u u 

^xy^yx 


(g 4 u € )6 + (G r )r 

v^y ^xy x^^y '^xy xy^ x; 


diL 


Combining terms in the equation above, obtain: 

- - - r E E u 

» = 1 — ^ 


xy'-y X 


X y '■^y X 


1 “ U U 

^xy^yx 




E u 

y ^x y 


1 " U. U 
^xy^yx 




+ G r 

*y xy 


di 


(2.40) 


The number of terms in equation 2.40 may be simplified using 
equation 2.39 to obtain: 


0 = ,in [-^ 

U L 1 - U 


X y^y X 


2 E u 

6 + 1 € € 

X y X 

1 “ u u 

'^x y y X 


1 " U U 

^xy^y: 


€: + 


G r 

X y X y 


dY. 


(2.41) 


The first variation of the total strain energy is taken with 
respect to the strains as extremal functions, noting first 
that in principle: 
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6 (a^ ) = 2 a(6a) 

5(aB) = a(6Q) + Q(6a) 


Therefore , 
6 




* y y * 


(2.42a) 

(2.42b) 


2 E,^yx 

^ ~ ^xy^yx 


€,6€ 


X 


2 E u ' 

x^yx 

1 - UL U 

^xy'^yx 


2 E, 


6 6e + 

X ^y 


1 ” U U 

^x y '^y X 


€,6£, 


+ 2 G r 6T 

xy xy xy 


dZ (2.43) 


The factor is cancelled in the process and again using 

equation 2.39 to change the third term in equation 2.43 to 

obtain the desired form for a recombination of terras; 

E E,«, 

- + 

1 - 

X y '^y X 


e'*>u = 


m 

V L 1 - 


y X 


1 " u u 

y X 




E u 

y ^x y 


1 - JJL U 

y ^y X 


e 6€ + 

X y 


1 ” u u 

'^x y ^y X 




+ G r 6 r 

xyxyxy 


(2.44) 


6“>U = 


HI 

V E I - n.,u,. 




1 “ Ul u 

^x y ^y X 


(^xy^x ^y^ ^^y 


+ G r 6F 

xyxyxy 


dY. (2.45) 
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Using equations 2.38a-c, the first variation in the total 
strain energy is expressed as: 

= III 'S -IS) 

V 

To express the first variation in the total strain 
energy in terras of stress and displaceraent , substitute the 
strain-displacement relations in equations 2 . 3a-c into equa- 
tion 2.46 to obtain: 


= JJJ I ^ ^ ^ [w%]2) 

V 

+ ay6(v0y - z wOyy + [w 0 ^] 2 ) 

+ T,y6(u0y + vO^ - 2 2 wO^y + wO^ wOy) ] di (2.47) 
Noting that the delta operator and differential operator are 
commutative as shown by the following in principle: 


6(a ) 

1 * 

= (6a). 

(2, 

, 48a) 

6 (a ) 

» 7 

= (6a), 

1 7 

(2. 

. 48b) 


Using equations 2.42a-b and 2.48a-b, equation 2.47 becomes 

■ JJJ C - 2 

V 


+ a {(6v0) - 2 (5w0) + wO (6w0) } 

J »/ *// >/ 1 / 


+ T^y { (6u0 ) y+ (6v0)^- 2 2 (5w0)jjy+ w%(6w°)y+ 

(2.49) 


Expressing the volume integral in equation 2.49 in terms of 

* 

plate thickness and the area of region Q , or middle surface 
of the plate (x-y plane), results in: 
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= II I - 2 (6wO)^^ + w?x^5wO)^} 

a -HH 

+ a {(6vO) - z (6wO) + wO (5wO) } 

/ * * 1/7 » / I 7 

+ T^y{(6uO)y+ (6vO)^- 2 z (6wO)^y} 

+ T__{wO (6wO),+ wO (6wO) } ] dz dx dy (2.50) 
The middle surface displacements u<> , vO , and w° are 
not functions of z. Therefore those displacements may be 
placed outside of the integral over the thickness of the 
plate. The result is 

♦ 5iH +>4H 

6“'a = Jj [{(6u0). 

o ->SH -HH 

♦ J4H +>4H 

+ {(6v0)y + w0y(6w0)y} J Oy dz - (6«°)yyj cTyZ dz 

-HH -HH 

+ H H + 5^2 H 

+ {(6u0)y+ (6v0)^> I T,y dz - 2 (6wO)^J T^yZ dz 

-HH -HH 

♦ HH 

+ {w0jj(6w°)y+ w®y(6w®)„} J T^ydz ] dx dy (2.51) 

-HH 

Substituting equations 2.14a-f, the force and moment result- 
ants, into equation 2.51 results in 

6“’o = n [ N.{(SuO). ♦ wo.(6v.O).} - M.(5wO)„ 

n 

+ Ny{(6v0)y + «0y(6w0)y} - My(6w0)yy 

+ N,y{(6u0)y+ (6v0)^} - 2 M,y(5w0)^y 

+ Njjy{w0^(6w0)y+ w0y(6w0)jj> ] dx dy (2.52) 
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The first variation of the potential of the applied 


forces and the normal load distribution is 


6^ ^ SuOjjds + ^ Nj^gSuOgds “ JJ ^ 6wo dx dy (2.53) 

T2 T2 q 

Note that N^^ is taken as positive in compression as shown 

in Figure 2.8. 

Using equations 2.52- and 2.53 to form the first varia- 
tion of the total potential energy it and setting that varia- 
tion equal to zero to satisfy the necessary requirement for 


extremization of tt, leads to 


6 7C 






6^'^ = IJ C - M*(6«°)xx 

Q 

+ N {(6v0) + wO (6w0) } - M (6w0) 

/ >7 >7 » 7 7 » 7 7 

+ N,,{(6u0)^+ (6v0)^} - 2 M,y(6w0)^^ 

+ Njjy{w0jj(6w0)y+ w®y(5w0)jj} ] dx dy 

+ ^ SuO^ds + ^ N^gSuOgds ~ J| ^ dx dy 

T2 T2 O 

= 0 (2.54) 

Referring to Figure 2.9, and using the chain rule for 
differentiation to establish a relation between the partial 
derivatives in x-y coordinates and the partial derivatives 
in n-s coordinates, results in 


^ lx 

= ^^Ix ( In 

^ ^=|x ( Is 

(2, 

. 55a ) 

^ |y 

= (^|y ( In 

^ ^=|y ^ 1- 

(2, 

,55b) 


Using Figure 2.9, equations 2.55a-b can be re-written to be 
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(2.56a) 


^ }j = ^ny( )n ]n (2.56b) 
From the geometry of Figure 2.9, considering n and s to be a 
rectangular set of coordinates at a point on the boundary 


a^jj = cos(n,x) = cos 0 = 
a^jj = cos(s,x) = cos(x/2 
a^y = cos(n,y) = cos' 0 = 
a.. = cos(s,y) = cos 9 = 

® 7 

where by definition: 
cos(n,x) is the cosine of the 
cos(s,x) is the cosine of the 
cos(n,y) is the cosine of the 
cos(s,y) is the cosine of the 
From equations 2.57a-d 


a = 

- a 

8 X 

ny 

a. « = 

a 


n X 


Using equations 2.58a-b, 
as 







(2.57a) 

+ 0) : 

= - sin 0 

- 

(s) 

f 

X 

(2.57b) 

sin 0 

II 




(2.57c) 

(s)^ 





(2 . 57d) 

angle 

between 

the 

n 

and 

x axes 

angle 

between 

the 

s 

and 

x axes 

angle 

between 

the 

n 

and 

y axes 

angle 

between 

the 

s 

and 

y axes 


(2 . 58a) 
(2.58b) 

2.56a-b can be re-written 


equations 



) n “ ®ny ^ ? s 

(2.59a) 

^ " ^ny^ 

}n + ]n 

(2.59b) 

the geometry 

of Figure 2.9 


a_, = cos 0 = 
n X 

dy 

3? 

(2.60a) 

a^jy = cos 0 = 

dx 

■ 3? 

(2.60b) 


In the x-y coordinate system, u<> and vO are the middle 


surface displacements. In the n-s coordinate system, 
and u®. are the middle surface displacements. From the 

B 


42 


uO z 

n 


uO 


^ny 

< 

o 

uO = 

s 

- 

uO 

+ 

^nx 

o 

> 


u<> = a_, uO„ - uO. 

n X n n 7 8 

vO = a__ u®_ + a__ u°_ 
n 7 n n X 8 


geometry of Figure 2.9, transformation from the x-y 
coordinate system to the n-s coordinate system yields 

(2.61a) 
(2.61b) 

Transformation from the n-s coordinate system to the x-y 
coordinate system yields 

(2.62a) 
(2.62b) 

Green's Theorem for integration by parts in two dimen- 
sions may be expressed in general as follows: 

J| a S^.dx dy = - JJ a^,S dx dy + ^ (aS) a^^^ds (2.63a) 

JJ a 0_ydx dy = - jj a yQ dx dy + ^ (o0) a^^ds (2.63b) 

Applying Green's Theorem as expressed by equations 2.63 
to change the terms in equation 2.54 which involve partial 
derivatives of the variations in uo , vO , and wO , and to per- 
mit extraction of those terms from the area integrals, re- 
sults in the following: 

J|^ Nj,(6u0)jjdx ~ dx dy + ^ Nj,6u0 

J| N,w0^(6w0)^dx dy 

= - J|^ (NjjWO JJ ) wO dx dy + ^ NjjW°jj6wO a^^j^ds 

1| .dx dy = - J| M.,.<6M0).dx dy M.(6v.o ),a„.ds 

dx dy - ^ M.,.®"” »n,ds + ^ <6»® ) ,a„,ds 

^2 ^2 

Note that Green's Theorem is applied for two successive in- 
tegrations to extract the second-order partial derivative of 
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the variation in wO whereas one integration is required to 
extract first-order partial derivatives of the variations in 
u° , vO , or wO . 

Ny(6v0)ydx dy = - y6v<> dx dy + ^ Ny6v0 a^^ds 

J| NyW0y(5w0)ydx dy 

= ■ n dx dy + ^ N,w0^6w0 a^^ds 

n r 2 

i| dy = - «,.,(6“'’>rdx dy +J M, (6 wO j 

= J1 “r.7T®”° dx dy ^ M, ,6w0 a„,ds + ^ (6 wO ) ^a„,ds 

“ ‘'2 f'2 

|| N.,(6u0),dx dy = - j| dx dy 

J| N^y(6v0)^dx dy = - N*y,,6v0 dx dy + ^ ^nx^^ 

2 

J| N,yW0^(6w0)ydx dy 

= - || (Nxy«?x^yS«'’ dx dy N^yWO^SwO a^^ds 

1| N,y«?y(6«°)xdx dy 

■ " J| ^^xy”?y^x®”° dx dy N^yWOy6wO a^^ds 

l[“xy(5«°^ydx = - I| «xy.x^5w0)ydx dy M,y(5w0)ya„^ds 

2 

= 11 dx dy - ^ M., ,6»o a„,ds ♦ ^ M,,, (5 wO ) ^a„,ds 

“ ^*2 ^2 

The equation immediately above was integrated first with 
respect to x and then integrated with respect to y. Revers- 
ing the order of integration provides a result which may be 
used additionally in view of the factor 2 in equation 2.54 

for the term involving the moment resultant M,,. 

* y 
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,y(6wO)^ydx dy = - 


II ”xy.y(««°>xdx dy + ^ M^y(6wO)^a„yds 
“ ^2 

= II ”x7.yx8«° dx dy - (f 6wO a„,ds + (f M ( 5 wO)^a ds 
^ t *2 1-2 

The results of having applied Green's Theorem to equation 
2.54 may be combined and written as follows for extreraiza- 
tion of the first variation of the total potential energy. 

" dx dy + ^ Njj6u0 a^^j^ds 

- dx dy + ^ ^x'^?x®”® ^nx*^^ 

- ^ ^y"? y ^ , y®'^® dx dy + ^ NyW0y6w0 a^y*^^ 

-II 


Ny.ySvO dx dy + ^ Ny6v0 a^^^ds 

^ 2 


^xy,y®u® dx dy + ^ Nj,y6u® a^^ds 

*' 2 


I| Njjy x^vO dx dy + |i N^y6v0 a^^j^ds 
I| ^^xy”?*\yS«® dx dy N,yW%Sw° a^yds 

I| ^^xy”?y) ,x5«® dx dy N.^wO^ewO a^^ds 

If Mx,xx««" dx dy + f »x,x««" ^nxds - § M,(6w0)^a„,ds 
^ ^2 ^2 

I| “y.yy®”® ‘'^!f ”y.y®"® ^ny^s My (6w0 ) ^a^^ds 

II Wxy.xy6«® dx dy + ^ M.y.xS”^^ a^yds - § M,y(6w0) a„,ds 

^ ''2 r -2 

II “xy.yxS”® dx dy + ^ M,y_y6w0 a„,ds - § ( 5 wO ) ^a^^ds 

« f 2 r*2 


i 


N_ 6 u 0 „ds + 

n n 


^ N„. 6 u«.d 3 - II 


q 6 w°dxdy =0 ( 2 . 64 ) 


a 
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Equation 2.64 is expressed in terms referred to the x-y 
coordinate system and provides a basis with which to deter- 
mine plate equilibrium equations and the associated boundary 
conditions. The equation is complicated but may be simpli- 
fied if a selection of terms is expressed with respect to 
the n-s coordinate system. 


Stress components in the x-y coordinate system trans- 
form as a second-order tensor to the n-s coordinate system 
as expressed by the following equation: 



{ na ) 

■ ^pj 

i 

( *y ) 

P.k 

=n,s j,i=x,y 

(2.65) 

The 

stress 

component 

"^nn 

is 




^nn 








n 

= ®nx 


^xx + 

^nx 







+ 

a a 

n y n X 

^yx ®ny ^ny ‘'"yy 

(2.66a) 

The 

stress 

component 

a 

is 




■Tna 

J 

^al 






5 

= ^nx 

^«x 

Txx 

X y 

^.y 



^ ^ny ®sx ^yx ^ ^ny ^ay ^yy f2.66b) 
The stress tensor is symmetric (equation 2.32) so that equa- 
tions 2.66a-b can be re-written as follows: 


n 

= a^ T 2 a 

n X XX X 

^ny 

^*y 

+ a^ 

ny 

^yy 

(2.67a) 


= a « a T > 

a 

a 

T 



n 8 

n X ax XX 

n X 

«y 

' *7 





^ny 

^ax 

T + 

X y 

a a T 

ny ay ' y y 

(2.67b) 


Using the following relations 




z a 

T = 

n n n 

X X 

X 

7 7 

= - a 

a 

= a 


n y ax 

n X 

sy 
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(2.68a) 


Equations 2.67a-b may be expressed as follows: 

ny y 



= a' a 

nx ^x 

+ 2 


^ny + 


= - 

^ny 


+ a^ T 
n X X y 





2 

” ^ny "*”xy 


+ ^ny ‘^y (2.68b) 

The force and moment resultants in the n-s coordinate 
system may be expressed in terms of stress by the following 
equations: 


Nn 

♦ HH 

dz 

(2.69a) 

Mn 

N„. 

+ >SH 

"-h1 

dz 

(2.69b) 

Mn. 

The 

force 

resultants in 

the n-s 


z dz 


na ~ J ^na 

->4H 


z dz 


(2.69c) 


(2.69d) 


expressed in terras of the force resultants referred to the 
x-y coordinate system by integrating equations 2.68a-b over 
the plate thickness and using the definitions given by equa- 
tions 2.69a-b and equations 2.14a-c which yields the follow- 
ing force resultant transformation equations: 



= a^ N + 2 a 

®nx X ®nx 

^ny 

«.y 

ny y 


(2.70a) 

Nn. 

= • a a N + 

nx ny x 

a^ 

n X 

«Ky 

- N., + a„, 

n y X y n x 

"ny ^ 

(2.70b) 

The 

moment resultants 

in 

the 

n-s coordinate 

system 

may be 


expressed in terras of the moment resultants referred to the 
x-y coordinate system by multiplying equations 2.68a-b by z, 
integrating over the plate thickness, and using the defini- 
tions given by equations 2.69c-d and equations 2.14d-f which 
yields the following moment resultant transformation equa- 
tions : 
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«, * 2 a„, a„, M,, + a^, M, (2.71a) 

Ma. = - a-. M. + a". ' a^, + a„. a„ M, (2.71b) 


= ®n. 6uO. 

n X n n 7 8 

6vO = a 6u° + a__ 6u<>, 

n 7 n n X 8 


Proceeding with a selection of terms in equation 2.64 
to be referred to the n-s coordinate system, and first con- 
sidering the boundary integrals, equations 2.62a-b may be 
used to express variations in uO and v<> as follows: 

(2.72a) 
(2.72b) 

Selecting boundary integrals which include variations in u° 
and vO from equation 2.64, those integrals are referred to 
the n-s coordinate system as follows: 

f N,6 u 0 a„,ds + ^ NxySviO a„yds 
**2 ^2 

+ ^ N^y6v0 a^^ds + ^ 6v0 a^^ds 


- f (N^a^, + N,ya„y)(a„,6u0^ - a„y6u0, 


) ds 


^ + Ny a„y)(a„y6u0^ + a„,6u0^) ds 



2 N,ya„a„, + N^aJ^) 6u0„ds 


^ 2 


+ if (- N. 
■■2 

^nx^ny ^xy^nx ~ ^xy^ny ^y ®nx^ny^ 

6u0^ds 

= ® NjjSuOjjds 

+ f N„,6u0 ds 

(2.73) 

1*2 

r 

‘ 2 



Equations 2.70a-b were used to obtain the final result. 

Equations 2.59a-b may be used to express derivatives of 
wO as follows: 


(wO) 


; X 


a ( wO ) - a ( wO ) 

n * ' ' n “ny ; 


(2.74a) 
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(wO)y = a„y(wO)^ + (2.74b) 
Selecting boundary integrals which include derivatives of 
and variations in w® from equation 2.64, those integrals are 
referred to the n-s coordinate system as follows: 


f N.w0^5w 0 a„^ds + § N,,yw0^6w0 a„yds 

+ f a„^ds + ^ NyWOySwO a„^s 

^2 ^2 

= § ^N,a„ + N,ya„y)(a„,(wO)^ - a„y(wO)^) 6w0 ds 

^2 

+ 6 (N a + N a )(a (wO) + a (wO) ) 6w0 ds 
j '”xy“ nx ‘y^ny'^ ny' (n **nx' 'a' 

^2 

= § (N.aJ. t 2 N.,a„,a„, * N,aJ, ) ( wO ) ^6W> ds 

I', 

di +N 3^ ~N 3^ Si )(w^)6w^ds 

X ' " x"nx“ny ^xy^ny ^y^nx^ny^'^ { ^ 

^2 

= ^ ds + ^ Nj^g(w0)g6w0 ds (2.75) 

^2 ^2 

Equations 2.70a-b were used to obtain the final result. 

Equations 2.62a-b may be used to express variations in 
wO as follows: 

(2.76a) 
(2.76b) 

Selecting boundary integrals which include derivatives of 
the moment resultants and variations in wO from equation 
2.64, those integrals are referred to the n-s coordinate 
system as follows: 

§ M,.x6w0 a^^ds + f Mxy,y5w<» a„,ds 


6w0 = a„. 6w0 

- a 6 w® 

nx n 

n y 

6w0 = a„, 6w0„ 
n y n 

+ a„, 6w0 
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^ My.y6w° a^^ds 
^2 ^2 

— V (M di + M a ) 6w^ ds 

J '“x,x®nx “xy,y®nx' 

'’2 

§ fMxy.x^ny ”y,y^ny^ ®wO ds 
^2 

= i [ (M,,, + a„, + a„y] 6wO ds 


Multiplication of the three-dimensional equilibrium 
equations of elasticity by z and then integrating over the 
plate thickness, in the absence of body forces, establishes 
relations among the quantities , M^ , M^^ , , and , 

where the latter two quantities are defined as the trans- 
verse shear force resultants from classical plate theory. 
The resulting equilibrium equations are: 

M*,x + M«y.y = Q, (2.78a) 

M*y.z My.y = Qy (2.78b) 

Qx.x Qy.y = 0 (2.78c) 

where q(x,y) is a normal load distribution on the plate. 

The right sides of equations 2.78a-b may be substituted 
into the boundary integral immediately above . 


C ^^x,x ^xy.y^ ^nx ^^xy . x ^y , y ^ ^ny] 

’2 

= f C Qx ^nx Qy ^ny ] ^S 

f '2 

= f [ Qn ] ds (2.79) 


Equation 2.61a, which transforms quantities in the x-y coor- 
dinate system to the n-s coordinate system, was used to 
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obtain the final result given in equation 2.79. 

Equations 2.59a-b may be used to express derivatives of 
variations in wO as follows: 


(6w0) 

> * 

= a (6w0) 


(2.80a) 

(6w0) 

» 7 

= a (5 ) 

+ a„^(6w0)^ 

(2.80b) 


Selecting boundary integrals which include the moment re- 
sultants and derivatives of variations in w® from equation 
2.64, those integrals are referred to the n-s coordinate 
system as follows: 


- f ^ M.,(6«0)^a ds 

^2 ^2 

- ^ M,jy(6w0)ya„^ds - ^ My(6w0)ya^yds 

3 - f (M,a„, + M,,a„,)(a„,(5»li;„ - a„,(6w0)^) ds 

^2 

- f <“.,a„. + «,a„,)(a„,(6»")„ + a„,( 6 « 0 )^) ds 

^ 2 

= - ^ <M.a = , t 2 * M,a^, H6»o )„ds 

^ 2 

- J <- «.a„,a„ + M„a^. - M.,a^, ♦ M,a„,a„, ) (6w0 ; ^ds 
^ 2 

= - ^ (6w0)^ ds - ^ (6 m°). ds (2.81) 


Equations 2.71a-b were used to obtain the final result. 

The second integral in equation 2.81 can be considered 
further. If the moment resultant is continuous and the 
boundary is smooth, then an integration by parts yields 

J Mn«(5w0)^ds = - J ,6w0 ds + [ M„,5w0 ] |J (2.82) 

1 X 
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For the case where the boundary is a smooth closed curve, 
the expression in brackets is equal to zero and the result 
obtained in equation 2.81 may be written as follows: 


- £ (6w0)^ ds - f (6w0)^ ds 
^•2 ’ ^2 • 

= - f (6w0)^ ds + f Mns.s^wo ds (2.83) 
^2 ^2 

For the case where the boundary is not smooth, such as a 
plate with four smooth sides which terminate at corners, or 
discontinuous points, then the expression in brackets in 
equation 2.82 must be evaluated over each section of smooth 
continuous boundary . 

Having simplified the boundary integrals, equation 2.64 
may be rewritten as follows for extremization of the first 
variation of the total potential energy. 


- [(«.“?. ).x^ (N., (N,w%),,]6«0 dx dy 

- Px.x* *'xT,,]®'^° [^XT.X* '',.,1®''° 

* ^ [M.,... * * N„. ♦ On] 6““ 

^ 2 

+ ^ (N„ + N„)6u0^ds + ^ + N^,)6u0^ds 

^2 ^2 

- <f M (6w0)^ ds - [ M„g6w0 ] ( =0 (2.84) 

where the last term must be evaluated over each portion of 
the boundary between discontinuities. If the boundary is 
smooth and continuous, that term equals zero. 
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Having satisfied the requirement for extreraization of 
the total potential energy by equating its first variation 
to zero, resulting in equation 2.84, a basis has been estab- 
lished for a series of deductions to be made. The middle 
surface displacements uO , v® , and w® , were specified to be 
arbitrary in the region Q. This means that variations in 
the surface displacements are not equal to zero in the re- 
gion Q. If equation 2.84 is to be satisfied, then those ex- 
pressions in brackets in the surface integrals which are 
multiplied by 6u®, 6v0, or 6w0, must be equal to zero. 
Those expressions are the equations of equilibrium and are 
written as 


N + N =0 

* , X X y , y 

N + N =0 
‘ X y , X " y , y 


M.,..* 2 (N.wo,) 

* * <J = 0 (2.85) 

( 2 . 86 ) 
(2.87) 

Equations 2.86 and 2.87 can be identified with the 
equilibrium equations for plane stress. Equation 2.85 may 
be simplified by using equations 2.86 and 2.87 to eliminate 
some terms when its product terms are expanded. Expansion 
of the product terms are written as 


(«x 

wO 

. X 

= 

N. w< 
* 1 

) + 

1 X X 


i , K 

«?X 

(2. 

. 88a ) 


wO 

> 

X ^ , y 

= ^y 

wO 
. xy 

+ 

Nxy 

.y «?x 

(2. 

,88b) 


wO 

1 

y>.x 

= ^y 

wO 
. yx 

+ 

«xy 

. X y 

(2. 

88c) 


wO 

. y 

^ » y 

Ny WC 

^ "f 

. y y 


f . y 

wO 

. y 

(2. 

,88d) 


Terms from each of equations 2.88a-b and equations 2.88c-d 
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can be combined and eliminated on the basis of equation 2.86 
and equation 2.87 as follows 

(N*,* ^ wO, = 0 (2.89) 

(Nxy.x + Ny.y) «?y = « (2.90) 

The terms in equations 2 . 88a-d which remain are 


N wO 

”x ",XX 


N wO 

*7 .*7 


N wO 

"xy ’^.yx 


N wO 
"7 ,77 


and the second and third terms may be combined. 
Equation 2.85 can be rewritten as 


M + 2 M + M 

“x , XX “xy . xy y , y y 


+ N.wo^. * 2 • q = 0 (2.91) 

This equation of equilibrium, opposed to the linear differ- 
ential equation of classical plate theory, is nonlinear. It 
contains the nonlinear terms consisting of products of the 
force resultants and curvatures. 

In addition to the means to determine the equations of 
equilibrium, equation 2.84 contains the natural (or force) 
and the kinematic (or geometric) boundary conditions which 
must be satisfied or specified on the boundary F 2 . These 
boundary conditions are 


Either 

n n 

or, u<Jjj 

is 

specified 

(2.92a) 

Either 

Nn. = -Nn, 

or, uo^ 

is 

specified 

(2.92b) 

Either 

Mn = 0 

or, wo^ 

is 

specified 

(2.92c) 

Either 

Q + M 

n a , a 





+ 

N„w0n+ N^^wO = 0 

n , n n a , a 

or , wO 

is 

specified 

(2.92d) 

And , at 

discontinuities [ 

M„.6w0 3 = 

0 


(2.92e) 


The last three conditions contain those terms which are 
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characteristic of conditions embodied in classical plate 
theory. However, the effective shear force condition in 
equation 2.92d, consisting of the shear force resultant term 
and the moment resultant term, has been made more complicat- 
ed by the addition of product terms (in-plane force result- 
ants multiplied by the slope at the boundary) . 

The equilibrium equations 2.06, 2.87, and 2.91 can be 
written using equations 2.25 and 2.26 as follows: 


Pii^S 

+ 

2^y 


(^I6:*2y],i * Pie^S 

+ 

Aye^y 

* *66rSy] 

1.7 = 0 









(2.93) 



^2 6^y 


A6«r»,].,+ Piz^S 


*22^? 

^26^Sy] 

= 0 
.7 









(2.94) 

Pi PS 

+ 

D 12 KO 

+ 

:^16*^Sy] , ,v 







2 Pie 

KO 

X 

+ d„k; * D,,K0,] , 

xy 





* ^ O33KO + Dj,K0^] „ 

* 2 [A,3€0 + 

+ [Ai 2«2 * *22^? + + q = 0 

(2.95) 

The equilibrium equations 2.93, 2.94, and 2.95 can be 

written in terms of the surface displacements using equa- 
tions 2.5a-c and 2.6a-c in substitution resulting in 
[A,,(U%+ + Aj3(v0^+ I4(w%)2} 

* IS(»%)2) * A26<''?,* 

* = 0 (2.96) 
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X 


[A^g{uO^+ }$(W0^)2} + A26(V%+ ^(W0y)2> 

■" A86<«%+ vO^+ «?x"?y>] , 

[Ai2{U%+ J^(wO^) 2} + A22(vOy+ >5(wOy)2> 

+ A2g{u0 + vO + wO wO }] , = 0 (2.97) 

CAll{u?X+ »5(w%)2} + A^2 {v%+ »5(wO^)2} 

+ A,g{uO^+ V%+ W?*W?y>] W?xx 

+ 2 [Ajg{uO^+ >$(w%)2} + A2g{vOy+ >5(wOy)2} 

+ Agg{uOy+ vO^+ wO^wOy}] wo^y 
[Ai 2^^?X-*- + A22(v0y+ J$(wOy)2} 

+ A 2 g{uOy+ vO^+ W?,W?y>] «?yy 

E®ll”?xx ■*■ ^X2”?yy ■*■ ^ ®ie”%y],xx 

2 [D^gWOjjjj + D26”?yy ^ ^6e”?xy],xy 

Pl2«?xx ^22^%y * 2 D 2 sWO^J 

+ q(x.y) = 0 (2.98) 

The laminate extentional stiffness properties A. j and 
the laminate bending stiffness properties D^j are functions 
of the spatial coordinates in equations 2.93, 2.94, and 
2.95. These equations are partial differential equations in 
X and y and may not be solved exactly. A finite element 
method will be used to obtain solutions for these equations. 
This will be discussed in detail in Chapter III. 
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CHAPTER III 


PROBLEM FORMULATION AND METHOD OF ANALYSIS 
FOR TENSILE RESPONSE 


Plate Geometries 

The geometric configuration of the plate under consid- 
eration is illustrated in’ Figure 3.1. The plate width is 
designated by W , the plate length by L , and the plate hole 
diameter by D . The positive horizontal x-axis is to the 
right and the positive vertical y-axis is 90° counterclock- 
wise from the x-axis . The plate is loaded by in-plane ten- 
sile forces in the x-direction distributed uniformly over 
each end. 

Four discrete plate geometries are considered. Two 
values of D/W equal to 1/3 and 1/6 are used along with two 
values of L/W equal to 1 and 2. The reason for considering 
various geometries is to ascertain whether an improvement in 
structural efficiency with the curvilinear fiber design is 
or is not a function of plate geometry. 

Equations Governing In-Plane Tensile Response 

A thin symmetrically laminated plate in a state of 
plane stress is considered. If the plate is subjected only 
to an in-plane tensile loading, the displacements developed 
due to that loading are only in-plane displacements. The 
out-of-plane displacements equal zero. For that situation, 
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Figure 3.1 Plate Geometry And Loading 
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equations 2.96 and 2.97 reduce to 


[Aiiuo^ + A,e(uO + yO )] 

+ Ase(uO^+ ^?«)].y 

[Ajguo^ + ^2gv0^+ Agg(uO + yO ^ 

^ A22^?y-*- A2 s(u%+ V%)] ^ 


= 0 


= 0 


(3.1) 

(3.2) 


These equations are lihear partial differential equations 
and the coefficients A^j (iij = 1,2,6) are functions of 
the coordinates x and y. These equations can not be solved 
exactly. A finite element method will be used to determine 
solutions. The type of boundary conditions considered are 
those associated with simply supported plate edges. 


Finite Element Model, Finite Element Meshes, And Boundary 
Constraints 

The plates are analyzed using a finite element code 
written previously to this study. As shown in Figure 3.2, 
the basic element is an eight-node isoparametric element 
with nine Gauss integration points. A relationship is need- 
ed between the element displacements at any point within the 
element and the element nodal point displacements , since re- 
sults are obtained with reference to the nodal point posi- 
tions. This is accomplished directly through the use of in- 
terpolation functions which are defined in a natural coor- 
dinate system. Coordinates of any point in the element in a 
global x-y coordinate system may be related to coordinates 
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Figure 3.2 Eight-Node Isoparametric Element 
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of the nodal points through functions expressed in a natural 
r-s coordinate system. 

The element shown in Figure 3.2 lies in the global x-y 
coordinate system. A transformation between positions of 
the nodal points and positions of any other points in the 
element is given by the relations 

X = $^(r,s) x^ •. (3.3a) 

y = ,|^ «.(r,s) y. (3.3b) 

where are the interpolation functions, x and y are the 

coordinates of any point of the element, and x^ and y. are 
the coordinates of the element nodes. 

Interpolation of the element displacements is accom- 
plished in the same manner as interpolation of the element 
coordinates as expressed by 

u = ^1^ $^(r,s) u. (3.4a) 

V = ,|^ 4^(r,s) V. (3.4b) 

where u and v are displacements of any point of the element, 
and u^ and are displacements of the nodes. The same in- 

terpolation functions are used to interpolate the element 
displacements and the element coordinates. This is the 
basis for the isoparametric finite element formulation. The 
interpolation functions for the eight-node element shown in 
Figure 3.2 are 

fj(r,s) = 1/4 (1 + r)(l + s)(-l + r + s) (3.5a) 

« 2 (r,s) = 1/2 (1 - r2)(l + s) (3.5b) 

$ 3 (r,s) = 1/4 (1 - r)(l + s)(-l - r + s) (3.5c) 
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$^(r,s) = 1/2 (1 - r)(l - s2) (3.5d) 
$ 5 (r,s) = 1/4 (1 - r)(l - s)(-l - r - s) (3.5e) 
«g(r,s) = 1/2 (1 - r2)(l - s) (3.5f) 
« 7 (r,s) = 1/4 (1 + r)(l - s)(-l + r - s) (3.5g) 
«g(r,s) = 1/2 (1 + r)(l - s2) (3.5h) 


The function $ and its subscript in equations 3 . 5a-h corres- 
pond to the sequence and position of nodal numbers in Figure 
3.2. (Reddy [3], pp. 242-254). 

The finite element meshes for the 4 geometries consid- 
ered in the tensile analysis are depicted in Figures 3.3 
through 3.6. The plate geometry, loading and material pro- 
perties are symmetric about the midspan of the plate in the 
X and y directions, so only one quarter of the plate is used 
in the analysis. Each quarter plate mesh consists of 192 
elements with 12 elements placed around the hole edge. 

The boundary constraints for the finite element mesh 
are depicted in Figure 3.7. The nodes at the edge of the 
hole are traction free. The displacements vO in the y-di- 
rection equal zero for the nodes on the edge which corres- 
ponds to D/2 < X < L/2 and y = 0, and that edge is traction 
free in the x-direction. The displacements u<> in the x-di- 
rection equal zero for the nodes on the edge which corres- 
ponds to X = 0 and D/2 < y < W/2 , and that edge is traction 
free in the y-direction. Tensile forces are prescribed in 
the x-direction at the nodes on the edge which corresponds 
to X = L/2 and 0 < y < W/2, and that edge is traction free 
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Figure 3.4 One-Quarter Finite Element 
For The Case L/W = 1 And D/W = 1/6 


64 



Figure 3.5 One-Quarter Finite Element 
For The Case L/W = 2 And D/W = 1/3 
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Figure 3.6 One-Quarter Finite Element 
For The Case L/W = 1 And D/W = 1/3 
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in the y-direction . The edge which corresponds to 0 < x < 
L/2 and y = W/2 is traction free in the x and y directions. 

Procedure For Determining Fiber Directions In Curvilinear 
Fiber Format Layers 

As stated previously, the basic idea behind the use of 
the curvilinear fiber forinat is to orient the fibers in the 
principal stress directions. In general, the principal 
stress directions depend on the stiffness properties. If 
fibers were introduced into a plate of isotropic material, 
for instance aluminum, forming in effect a layer, its stiff- 
ness properties would be changed. However, its stress state 
in the thickness direction would not be changed necessarily 
and a principal stress direction could be determined. In 
contrast, determining principal stress directions for a 
laminate is not a meaningful or comparable concept. The 
response of an entire laminate is dependent on the stiffness 
properties of each layer which may be different layer by 
layer. The stress state in each layer of a laminate is 
different. Thus, principal stress directions may be deter- 
mined for each layer and those directions will generally 
vary from layer to layer. An iterative process was used to 
determine the fiber directions within a group of layers 
(designated curvilinear layers) such that the fibers were 
aligned everywhere with the principal stress directions in 
those curvilinear layers. Other layers (designated to have 
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a straightline fiber format at a specified angle), which may 
be included in a laminate design, would not be subject to 
any required realignment with principal stress directions. 
The iteration process determines the principal stress direc- 
tions only for the curvilinear layers. With the procedure, 
the fiber directions are specified relative to the positive 
horizontal x-axis. 

Due to use of the finite element method, two approxima- 
tions are inherent to the analysis. First, the fiber direc- 
tions are assumed to be constant within an element. Second, 
calculations for principal stress directions , and hence 
fiber directions, are keyed on the stress computed at just 
one Gauss point in the element. Gauss point number 8 was 
chosen in order to use stresses which are as close to the 
hole edge as possible. Figure 3.8 shows the location of 
Gauss point 8 as the particular eight-node element was used 
in this study. Gauss point 8 was chosen rather than Gauss 
points 7 or 9 because Gauss point 8 is the central point in 
that group of three and provides a better average value for 
the stresses. 

The iteration process is described in the following as 
a series of numbered steps in procedure: 

1. The iteration process may be initiated with the 
assumption that the laminate is an isotropic plate of alumi- 
num. Thus, only two independent material properties among 
the modulus of elasticity E, the shear modulu.s G, and the 
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The isotropic material property values may be used as the 
required input to the finite element program to provide 
results for laminate strains €y, and and laminate 

force resultants N^, N^, and at each of the 9 Gauss 

points in each element of the finite element mesh. As men- 
tioned previously, only the results associated with Gauss 
point 8 are used to establish the fiber direction character 
of each element. 

2. The iteration process continues with the assump- 
tion that the laminate now consists of layers of orthotropic 
material. Equations 2.9a-d are applied to determine values 
of the reduced stiffnesses Q, obtained with respect to the 
1-2 coordinate system. These values of reduced stiffness 
are the same for every element in the finite element mesh 
(as functions only of material properties). 

3. At this point, if the fiber directions for the 
curvilinear layers in the laminate were known, values of the 
reduced stiffnesses Q, with respect to the x-y coordinate 
system, could be computed for each element in the finite 
element mesh using equation 2.12. These values of reduced 
stiffness would be different from element to element of the 
finite element mesh. Initially, however, it is assumed that 
the fiber directions in all of the elements are equal to 
zero. Values of reduced stiffness with respect to the x-y 
coordinate system may be computed on that basis (zero direc- 
tion) . Recall that a laminate may consist of straightline 


71 



fiber format layers as well as curvilinear fiber format 
layers. The stresses for the curvilinear fiber format 
layers in the x-y coordinate system should be computed sep- 
arately from those of the straightline fiber format layers 
in order to determine a new set of element principal stress 
directions for the curvilinear fiber format layers. Values 

f 

of stress computed for one curvilinear fiber layer would be 
the same for every layer in a group of curvilinear fiber 
format layers within a laminate in the x-y coordinate sys- 
tem. It follows that the element principal stress direc- 
tions would be the same for each curvilinear layer within a 
laminate . 


4. Having values of laminate strain from step 1. for 
Gauss point 8 of each element and all fiber layers , and from 
step 3. , values of the reduced stiffness for the curvilinear 
fiber format layers, equation 2.13 is applied to compute the 
stresses in each element for the curvilinear fiber format 
layers within the laminate. 

5. The principal stress direction for the curvilinear 
layers in each element is computed using the equation 


0 = 1/2 tan*i 




(3.6) 


where the principal stress direction for an element is spec- 


ified by 9, and . Cy , and are the stresses computed 

in step 4. At this stage in the iteration process, it is 


convenient to designate these principal stress directions as 
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the current fiber directions for the curvilinear fiber for- 
mat layers. 

6. Applying equation 2.12, compute with respect to 
the x-y coordinate system, a set of reduced stiffness val- 
ues Q in each element for the curvilinear fiber format 

layers, using the values of reduced stiffness Q from step 2. 

* 

and the principal stress direction values obtained in step 
5. 

7. Applying equation 2,12, compute with respect to 
the x-y coordinate system, a set of reduced stiffness val- 
ues Q in each element for the remaining straightline fiber 
format layers, using the values of reduced stiffness Q from 
step 2. The fiber direction is specified and characteristic 
of each layer. For example, the fiber direction of a pair 
of layers may be specified to have fibers oriented at +45° 
and -450 relative to the x-axis of the laminate. 

8. Applying equation 2.20 to all of the layers in the 
laminate, compute a set of extensional stiffness values A^j 
for each element of the finite element mesh using values of 
the reduced stiffness obtained in step 6. and step 7. 

9. Divide the extensional stiffness values A. . , by 

^ J 

the value of the thickness of the laminate to obtain lami- 
nate stiffness property values for each element of the fin- 
ite element mesh. 

10. Input the laminate stiffness properties for each 
finite element into the finite element program. The finite 
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element program produces a new set of results for laminate 
strains Gy, and Tj^y, at each of the 9 Gauss points for 

each element. 

11. Using the values designated as current fiber di- 
rections for the curvilinear fiber format layers obtained in 
s'tcp 5. and values of the reduced stiffness Q for orthotrop- 
ic material obtained in step 2,, apply equation 2.12 to 
compute new values for the reduced stiffness Q . 

12. Using the values of reduced stiffness Q, obtained 

in step 11. and the values of laminate strain for Gauss 
point 8 obtained in step 10. , apply equation 2.13 to compute 
values for the stresses and r^y, in each element for 

the curvilinear fiber format layers. 

13. Using the values of stress obtained in step 12., 
apply equation 3.6 to determine the principal stress direc- 
tions in each finite element for the curvilinear layers. At 
this stage in the procedure, these principal stress direc- 
tions are now designated as current fiber directions 0c, for 
the curvilinear fiber format layers. The values for the fi- 
ber directions obtained in step 5. are now redesignated as 
previous fiber directions 0p. With this step, one cycle of 
the iteration procedure is completed. 

The fibers have been aligned with new (current) prin- 
cipal stress directions. The iteration procedure is simple, 
straightforward, and repetitive, and may be continued until 
the principal stress directions, and thus the fiber direc- 
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tions obtained previously, can be compared to current fiber 
directions within a chosen tolerance. That tolerance level 
is considered to have converged when the following quantity 
(relative error in 6c) is determined to be less than 0.01 
for each element. 

6 p - 0 c 

< 0.01 , ( 3 . 7 ) 

0P 

where 0p specifies the fiber direction for an element as a 
result of a completed previous iteration and 0c specifies 
the fiber direction for an element as a result of the com- 
pleted current iteration. 

Additional cycles in the iteration procedure may be 
accomplished in the following steps. 

1. Using values of the current fiber directions 0c 
for the curvilinear layers in each finite element. and val- 
ues of the reduced stiffness Q from application of equations 
2.9a-d, apply equation 2.12 to compute values of the reduced 
stiffness Q. 

2. Using values of the reduced stiffness Q for the 
straightline fiber format layers which have specified fiber 
directions, apply equation 2.12 to compute values of the 
reduced stiffness Q. 

3. Using values of the reduced stiffness Q obtained 
in step 1. and step 2. above, apply equation 2.20 to compute 
values of the extensional stiffness A. ^ for each finite ele- 
ment of the laminate. 
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4 . 


Divide the extensional stiffness values A . by the 
thickness of the laminate to obtain values of the laminate 
stiffness properties for each finite element in the finite 
element mesh. 

5. Input the stiffness properties obtained in step 4. 
into the finite element program to obtain values of laminate 
strain and r__. 

6. Apply equation 2.13 to obtain values of stress 

Oy , and T^y, using the values of reduced stiffness Q obtain- 
ed in step 1 . and laminate strains obtained in step 5 . 

7. Apply equation 3.6 to obtain the principal stress 
directions using the stresses obtained in step 6. 

8. Apply equation 3.7 to check if convergence of 
fiber direction values has been achieved within the chosen 
level of tolerance of relative error. 

The circumferential stress around the hole is tensile 
at the net-section and compressive at the horizontal center- 
line of the plate. The circumferential stress changes from 
tension to compression around the hole edge when proceeding 
from the net-section to the horizontal centerline of the 
plate. There is a point on the circumference of the hole 
where the circumferential stress equals zero. The hole edge 
is traction free so that all of the stresses equal zero at 
that point. That point is called the isotropic point at the 
hole edge. All directions are principal stress directions 
at the isotropic point so the fiber direction at that point 
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was chosen to be consistent with the fiber directions of 
adjacent elements. 

Elastic Material Properties 

The elastic material properties for a lamina considered 
in this study were 

Ej = 19.885 X 106 psi Ej = 1.281 x 106 psi 

Nu ^2 = 0.298 G ^2 = 1.003 x 1Q6 psi 

and the lamina thickness was chosen to be 0.005 in. 

These elastic material properties were chosen to repre- 
sent closely the commercial material AS4/3501. Relations 
from micromechanics based on the rule-of-mixtures were used 
to determine the elastic material properties. The relations 
are given by 


Ex = 



(3. 

8a) 

Ea = 

(1 + V* )/(1/E,2 

* ■'t /E») 

(3. 

8b) 

Nu,2 

= VjNuj + v^Nu^ 


(3. 

8c ) 

Gx2 

= (1 + v*2)/(1/G, 


(3. 

8d) 

: the 

parameters in equations 3 . 9a-d are defined 

as : 


< 

11 

0.65 

fiber volume fraction 



= 

0.35 

matrix volume fraction 



Nu, = 

0.27 

fiber Poisson's ratio 



= 

nt 

0.35 

matrix Poisson's ratio 



II 

> 

0.269 

reduced matrix/fiber volume 
ratio for E 2 


II 

CM 

> 

0.162 

reduced matrix/fiber volume 
ratio for G ^2 
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Eji = 30.32 X 105 psi fiber elastic modulus in fiber 

direction 

E ^2 ” 2.248 X 105 psi fiber elastic modulus perpen- 

dicular to fibers in the plane 
of the lamina 

E^ = 493 X 103 psi matrix elastic modulus 

Gj = 3.683 X 105 psi fiber elastic shear modulus 

= 182.59 X 103 psi matrix elastic shear modulus 

See Tsai [4], pages 11-2 thru 11-11 for a detailed discus- 
sion of these relations. 

Failure Strength Theories 

Two failure strength theories were used in this study 
to predict the failure load in a laminate. The first is a 
noninteracting strain-based criterion called the maximum 
strain failure criterion and the second is an interacting 
stress-based criterion called the Tsai-Wu tensor polynomial 
criterion . 

Maximum Strain Criterion 

According to the maximum strain criterion, five para- 
meters are used to determine how failure occurs in a lami- 
nate. These parameters are the failure strains for the ma- 
terial used. The notation is as follows: 

: Tensile failure in fiber direction 

: Compression failure in fiber direction 

€|.j. : Tensile failure perpendicular to fiber 

: Compression failure perpendicular to fiber 
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rjj • In-plane shear failure 

The five strains constitute five failure modes and are 
assumed to be independent from one another. The strength 
ratio R is determined by dividing the failure strain by the 
actual strain due to the applied load. There are five 


;trength 

ratios, one 

for 

each 

mode 

of 

failure . 


Rit 

= 


= > 

Used 

when 

^1 

> 

0 

(3.9a) 

®ic 

- 


= > 

Used 

when 


< 

0 

(3.9b) 

^2T 



= > 

Used 

when 

^2 

> 

0 

(3.9c) 

^2C 

= 

1 ^ i C 2 1 

= > 

Used 

when 

^2 

< 

0 

(3.9d) 

^12 

= 

1 2 /^ 1 2 1 







(3.9e) 


where Rj.j, is the tensile strength ratio in the fiber direc- 
tion, Rj(, is the compressive strength ratio in the fiber di- 
rection, Rj-j. is the tensile strength ratio perpendicular to 
the fibers, compressive strength ratio perpendic- 

ular to the fibers, and R ^2 is the shear strength ratio. 

For a given level of applied load, the strength ratio R 
which has the lowest magnitude among the five possibilities 
indicates the mode of failure for the layer in the laminate. 

Tsai-Wu Criterion 

The Tsai-Wu criterion takes into account that modes of 
failure can be coupled and interact with one another. There 
are six strength parameters associated with the Tsai-Wu cri- 
terion. These are F^^. F^, F^ . F^^. and F*^. F*^ is 

treated as an empirical constant and has a normalized value 


79 



of (-0.5). 

The strength to stress ratio R is determined from the 
quadratic equation: 

-b ± (b2 + 4a)>4 

R = (3.10) 

2a 

The values of a and b are determined from: 

a = F„af - 2F;,(F,-.F„)!i<,,aj + F„a| * F„Tf, (3.11a) 
b = F,o, + FyUj (3.11b) 

The stress values are stresses due to the applied load and 
the strength parameters depend upon the material used. 

Table 3.1 lists the failure strains for the maximum 
strain criterion and the failure stresses for the Tsai-Wu 
criterion . 


Table 3.1 Failure Strains And Failure Stresses 


Failure Mode 

Strain 

Stress 


(10-3) 

(103 psi ) 

Tensile failure in fiber direction 

10.49 

209.82 

Compression failure in fiber direction 

10 . 49 

209.82 

Tensile failure perpendicular to fiber 

5.77 

7.54 

Compression failure perpendicular to fiber 

22.99 

29.87 

In-Plane shear failure 

13.10 

13.49 


Application of the two failure strength theories to 
predict the failure load and failure mode of a thin laminat- 
ed fiber-reinforced plate will be discussed in Chapter V. 
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CHAPTER IV 


PROBLEM FORMULATION AND METHOD OF ANALYSIS 
FOR BUCKLING RESPONSE 

Equations Governing Buckling Response 

The equations which govern buckling of a thin flat 
plate subjected to in-plane loading only (q(x,y) = 0) are 

derived in this chapter. The method used to derive the 
buckling equations is the adjacent-equilibrium criterion in 


which the middle 

surface displacements are given by 


uo 

= u° + 

6u^ 

(4.1a) 

vO 

= vO + 

€v0 

(4.1b) 

wO 

= w° + 

€w0 

(4.1c) 


where (u°,vO,wO) and (u0,v0,w0) are equilibrium configura- 
tions, € is a small parameter, and (u^,v^,w^) are arbitrari- 
ly small incremental displacements. The thin flat plate 
considered here is subjected to in-plane edge loading only 
so wO and its derivatives equal zero. (Brush and Alraroth 
[5]. PP • 90-91). With that condition, equations 4.1a-c 
become 


uO 

= uO + €u^ 

(4, 

, 2a) 

vO 

= vO + €v^ 

(4, 

, 2b) 

wO 

= €w^ 

(4. 

2c) 


The equilibrium equations were derived in Chapter II. These 
equations, for the case where the normal load distribution 
q(x,y) equals zero, are 
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+ A,5{u0y+ V0^+ «?X«?y>] 

CAi6^'^?x'^ J$(wO^) 2} + A2s{v0y+ J5(wOy)2} 

+ A8g{u0y^• vO^+ ,7=0 

[Aj,{uO^-h >$(w0,)2} + A25{vO^+ >S(wO^)2} 

+ Ag,{uO^+ vO^+ wO^wO^}] , 
CAi2^'^?x'^ ^(w0^)2> + A22(v%+ >$(wO^)2} 

+ A2e(u%+ vO^+ wO^wO^}] y = 0 

*^(w%)2} + A^2(vOy+ >$(w0y)2} 

+ A^e{uOy+ vO^+ «?*«?y>] w?xx 

+ 2 [Aie{uO^+ J$(wOj2} + A2g{v0y+ J$(w0y)2} 

+ A8g{u0y+ vO^+ wO^wOy}] 

^ [Ai2{u0^-K ^,(w0^)2> H- A23(v0^+ V,(wO^)2} 

+ A 26 (u%+ vO^+ «?,«?y>] «%y 

[Oll”%x Oi2^°yy ■*■ 2 Dj gWO jjy] 

2 pj^gwOjjjj + D2gW0yy ■»■ 2 DggwOjjy] jjy 

- [Ol2«?xx 022 «%y + 2 D 2 gW 0 ^J_y^ = 0 


Substituting equations 4.2a-c into equations 4.3, 4.4 
4.5 results in 


6u0 + A,2(v0y+ Gvo y + 

+ A,g{ao^+ €uO y+ vO^+ €vO 62w0^^w0 

"■ [Ai6{^?x+ ’ 5 (€wO J2} + A2g{v0y+ GvO y 

+ Agg{Go^+ ^^?.y+ ^?x-" ^"”?,x«?,y>] . 


»^2(GwO y 

,y>].« 

+ »^(6wO 


(4.3) 

(4.4) 


(4.5) 
and 

2 } 

) 2 } 
y 

(4.6) 
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[A,g{ao^+ 6 u? + A 2 g{vOy+ >$(€w?_y) 2 } 

+ Agg{aOy+ eu 5 y+ vO^^- ev^_^+ € 2 «?,x«?,y>] .X 

^ [Al2(ii?x^ ,,)2} H- A22 {v0^+ €v? ^2(6w? y)2} 

A2e{uOy+ y+ vO^+ evj ,,+ €^«?.x«?,y>] ., = 0 (4.7) 

[A,,{Qo^+ €u? ,+ _J2} + A^2(v%+ €v} ,+ >^.(6W?^^)2} 

+ A,g{GOy+ eu? y+ vO^+ e^?.x-" ^ 2 «?,x«?.y>]^W?.xx 
+ 2[A,g{Go^+ 6u? J$(GwO J2} + A2e(v0^+ GvJ >$(6w0_y)2} 

+ Agg{GOy+ Gu?_y+ vO^+ ^2”?.x«?.y>]^«?.xy 

+ [A,2{G%+ eu?_,+ »$(Gw? _,)2} + A22 (vO^+ »^(6wO _y)2} 

+ A26{Go^+ ^^?,y-" ^?x^ ^^?,x-" ^^«?,x«?.y>]^”?.yy 
[|D^ j6wO xx"*" (^1 2^”? , y y‘‘‘ ^ 1 6^” I , x y ] , x x 

- D26Gw5^yy+ 2 Dg^ewO 

Pl2^”l,xx‘'’ (^2 2^”? , y y '*' ^ I , x y ] , y y ~ ^ (4.8) 

Grouping equations 4.6, 4.7, and 4.8, by powers of 6 results 


in 


[Aii^?x ^ A, 2 V% + A,g(GOy + 

+ ^ [AiiU?.*+ Ai 2 V?,y+ A^g(u? y + v^ 

+ G2[A,,»$(wO ,)2 + A,2’^(«?,y)^ A,g(wO 

CA XG^^x "*” A^b^^w Arb(^®' '*’ 


‘2G’':y " "GG^'^Ty " '^?x^] ,y 


X 

A 2 G 

V? 

, y"*" Age 

(^?.y 


, y 

?. 

x)^ 

+ 

A 26 *$(w? 

.y>^ 

* Aesf"?, 

• "?.y>] .y 

2 6 

^?y 

+ 

AggCuOy 

+ vO 
1 



+ 

X 

A 2 8 

vO 

, y"^ Age 

(^?.y 


, X 


x)^ 

+ 

A26^(«? 



."?,y>] ,. 
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[Ai2^?x AjjVOy + A26 (u% + ^?x)].y 

e [A 12 U? „+ A22V?,y+ A26(uO_y + J] y 

+ €2[A,2>^(w^ ,)2 + A22>^(wf^y)2 + A^ « ( W? ^ ^ ^ ^ = 0 (4.10) 

e [AijGo^ + A^jvOy + A^gCaOy -*- 

•" ^^[AiiU?,*-*- Ai2V?.y+ A,g(uJ y + 

^ e3[A,,J$(w? x)2 + A,2>^(wf^y)2 + A,g(w^ 

■*• 2G [A^gUO^ + A^gVOy + AggCGOy + ^?x)]«?,xy 
+ 262[A,gU?^,+ A2gVf_^+ Agg(u?^y + v?.x)]"?.xy 
-K 263[A,g>^(w? ,)2 + A2 ,»$(w^ ^)2 + AggCw? ^)]w? 

+ € [Ai2^?k ■" A22V% + A2g(u% + vojjwj yy 

"■ €''[Ai2^?.x-^ A22V?^y+ A2g(uO y + V^^,)]wO yy 

+ €3[A,2’^(«?.*)2 A22»^(w?.y)2 ^ A2g(w^ 

G [Di^wO ^12”?, yy'*' 2 Dis”i,xy],xx 

26 [Die^f.xx'*' ^26"?. yy"^ 2 DggWO 

- ^ Pl2«?,xx-^ ^22«?.yy^ 2 ^ ^ = 0 (4.11) 


Referring to equations 4.9, 4.10, and 4.11, the pre-buckling 
force resultants for the equilibrium configuration are 


Nx = 

AiiU% 

+ 


+ 

Aie(^?y 

+ 

vO ) 
} * 

(4. 12a) 

^ = 

Ai2^?x 

+ 


> 

A 26 <^?y 

+ 

^?x) 

(4.12b) 

= 

Aie^?x 


*26''?, 

+ 

Asg(il0y 

+ 

vO ) 
1 * 

(4 . 12c) 


Referring again to equations 4.9, 4.10, and 4.11, the buck- 
ling force and moment resultants are 


X 1 ~ 

An'J?,x‘^ 

Al 2 ^?.y^ 

Ais(^?,y 

+ 

^?.x) 

(4. 13a) 

y 1 ■ 

Ai2'^?.x + 

A 22 '^?,y'" 

A 26 , y 

+ 

V?,x) 

(4.13b) 

xy 1 = 

Aig^?.x-^ 

A26'^?,y-" 

Aes^^? , y 


V?.x) 

(4.13c) 
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(4.14a) 


~ “ ^12”?, yy'*’ ^ ® 16 ”l,xy] 

^yl " ■ Pl2”?,xx‘*‘ ^ 22 ”?, yy"^ ^ ^26”?, xy] 

^xyl ~ ~ Pie'^^.xx^ ^2s'^?,yy^ ^ ^Se'^l.xy] 


(4.14b) 

(4.14c) 


The terms which involve and vO in equations 4 . 9 and 
4.10 sum to zero in view of equations 2.86 and 2.87, equa- 
tions 2.96 and 2.97, and equations 3.1 and 3.2, all of which 
refer to the equilibrium configuration. Therefore, those 
terms may be eliminated from equations 4.9 and 4.10. The 
terms which involve the incremental displacements in equa- 
tions 4.9, 4.10, and 4.11 contain terms which are linear, 

quadratic, and cubic in the small parameter 6. Since £ is 
small, terras associated with quadratic and cubic values of G 
may be neglected. The terms which remain comprise the buck- 
ling equations (having omitted the small parameter €). 

[AiiU?,x+ A^s(u0 y + _,)] ^ 

^ [Ai6'"?.x-^ A2evj_y+ Ag3(u0 y + vO = 0 (4.15) 


[Ai6^?.x+ A^gv^ y+ Aeg(u0 y + v^_J] 

+ [Ai2^?,x-" A 22 V? y+ A26 (u? y + vO ,)] y = ( 

[Aji^o^ + A^2V% + A^e(G0y + vOj]„o^^^ 

2 [A,6U% + A2gv% + Aeg(aOy + vO^)]wO ,y 

[Ax2^?x ^ A22v0y + A2e(G0y + vojjwo yy 

Pll”l,xx'*' ^ 12 ”?, yy"*" 2 DjgWO jjy] 

- 2 Pi6'^?.xx+ I^26«?.yy"' 2 D6eW?.xy],xy 

Pl2”?,Xx‘‘" ^ 22 ”?, yy'*’ 2 D26”l,*y],yy 
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(4.16) 


= 0 (4.17) 


C-'A 



Substituting equations 4.12 to 4.14 into equations 4.15 to 
4 . 17 results in 


Nxyi.x-N,!., = 0 

^xl.Kx 2 ^xyl.xy ^yl.yy 


(4.18) 

(4.19) 


+ „ + 2 = 0 ( 4 . 20 ) 
Equation 4.20 is uncoupled from equations 4.18 and 4.19. 
Equations 4.18 and 4.19 are independent of the pre-buckling 
state or equilibrium configuration and are partial differen- 
tial equations. Their only solutions are u^ = 0 and vj = 0. 
Equation 4.20 is the governing partial differential equation 
for buckling analysis and in terras of displacements has been 
written previously as equation 4.17. Values of , when de- 
termined, provide a measure of the amount of buckling. The 
laminate extentional stiffness properties and the lami- 
nate bending stiffness properties ^ are functions of posi- 
tion in X and y in the plane of the plate. Equation 4.17 
may not be solved exactly. A finite element code using the 
Engineering Analysis Language ( EAL) [6] was used to obtain 
approximate solutions and numerical results. 


Finite Element Model Used In The Engineering Analysis 
Language (EAL) 

The finite element model used in EAL is formulated us- 
ing the assumed stress field - minimum complementary energy 
method developed by T, H. H. Plan [7]. The element used is 
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a four-node combined plate membrane and plate bending ele- 
ment having five degrees of freedom per node as shown in 
Figure 4.1. 

The strain energy of an element for linear elastic 
behavior in terms of the generalized displacements {q} is 

U = 1/2 |_qj[ k ]{q} (4.21) 

where [k] is the element stiffness matrix. 

The minimum complementary energy is 

= U - JJ = minimum (4.22) 

where U is the strain energy in terms of the stress compon- 
ents, u^ are the prescribed boundary displacements, and 
are components of forces on the boundary. These force com- 
ponents are related to stresses directly next to the boun- 
dary through Cauchy's formula 




(4.23) 


where n. is the direction cosine of the boundary normal. 

J 

The stress distribution is expressed in terms of unde- 
termined stress coefficients {Q} such that 

<^22 ^33 J 

{fl}^ = 02 03 ••• J 

{c} = [ P ] {Q> (4.24) 

The elements of the matrix [ P ] are functions of position. 
The stress-strain relations are 

{€} = C N ] {a} (4.25) 

Using equation 4.25, the internal strain energy is 

(4.26) 


0 = 1/2 11^ [aj [ N ] {a} d'£ 
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Figure 4.1 Four-Node Quadrilateral 
Combined Membrane And Bending Element 
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Using equation 4.24, equation 4.26 may be expressed as 

U = 1/2 |_SJ [ H ] {0} (4.27) 

where the matrix [ H ] is 

[ H ] = [ P ]^[ N ] [ P ] (4.28) 

The prescribed boundary displacements are 

{u} = [ L ] {q} (4.29) 

f 

where the terras in the (natrix [ L ] are functions of posi- 
tion on the boundary and the terms in the vector {q> are the 
generalized displacements at the nodes. Since the forces 
{S} on the boundary may be expressed in terms of the 
stresses {c} as in equation 4.23, these forces may be relat- 
ed to the undetermined stress coefficients {6} through 
equation 4.24 such that 

{S} = [ R ] {0} (4.30) 

where the terms in the matrix [ R ] are functions of posi- 
tion on the boundary. 

The total complementary energy given in equation 4.22 
may be written as 

= 1/2 l_0J [ H ] {0} - 1_0J [ T ] {q} (4.31) 

where: [ T ] = Jj [ R ] [ L ] dA (4.32) 

The condition to be satisfied in order for the total 
complementary energy to be a minimum is that the partial 
derivative of the complementary energy with respect to each 
of the undetermined stress coefficients Oi, (i = 1,2,3,... ) 

should be equal to zero. Using equation 4.31 and taking 
partial derivatives produces the following result 
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t H 3 {S} 


C T 3 {q} 


(4.33) 


The undetermined stress coefficients {0} may be obtained by 
pre-multiplying equation 4.33 by the inverse of [ H ] which 
produces the result 

{(3} = [ H ]-i[ T 3 {q} (4.34) 
Substituting equation 4.34 into equation 4.27 results in 

0 = 1/2 {q}T [ T 1T[ H ]-i[ T ] {q} (4.35) 
Comparing equation 4.35 to equation 4.21 provides a relation 
with which to determine the stiffness matrix [ k ] . 

[ k 3 = C T 3T[ H 3-K T 3 (4.36) 
Increasing the number of assumed stress coefficients in 
effect decreases the excess strain energy and provides a 
more accurate solution. Increasing the number of elements, 
which increases the size and number of stiffness matrices, 
also leads to an increase in accuracy of the results. 

Finite Element Meshes And Boundary Conditions 

The finite element meshes for the four geometries con- 
sidered in the buckling analysis are depicted in Figures 4.2 
through 4.5. Each upper right quadrant in these meshes is 
identical to the one-quarter plate meshes considered in 
Chapter III. In other words, the mesh configurations are 
the same in each chapter. There are 192 elements in each 
quadrant with 12 elements placed around the edge of the hole 
or 768 elements for the entire plate with 48 elements placed 
around the entire edge of the hole. The complete plate con- 
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Figure 4.2 Finite Element Mesh For Case 
L/W = 2, D/W = 1/6 
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Figure 4.4 Finite Element Mesh For Case 
L/W = 2. D/W = 1/3 

93 







figuration of 768 elements is used in the buckling analysis 
instead of only a one-quarter plate configuration because 
only one buckling analysis per design is required. Four 
buckling analyses per design would be required with the one- 
quarter plate configuration (one analysis per quadrant). 
Four sets of boundary conditions would be required in a one- 
quarter plate analysis ih order to obtain results for all of 
the possible buckling modes. The four sets of boundary con- 
ditions for the one-quarter plate configuration are symmet- 
ric-symmetric, symmetric-antisymmetric, antisymmetric-sym- 
metric, and antisymmetric-antisymmetric . The choice was 
made to use one set of boundary conditions applicable to the 
complete plate in order to obtain results for all of the 
possible buckling modes. The boundary conditions imposed 
are those for a simply supported plate. 

The following two chapters contain the results of the ten- 
sile and buckling analyses. 
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CHAPTER V 


TENSILE STRENGTH OF THE PLATE DESIGNS 
Numerical Results 

This chapter presents the results for the uniaxial ten- 
sile loading problem. The failure loads and failure modes 
of four plate geometries are considered. Each geometry is 
used in conjunction with four straightline composite layups 
and five curvilinear composite layups to form sixteen-layer 
laminate designs. Tables 5.1 through 5.4 show the predicted 
tensile failure loads for the 36 designs considered. These 
tables are arranged according to the plate geometries. The 
failure loads shown in Tables 5.1 through 5.4 are based on 
the maximum strain criterion and the Tsai-Wu criterion. The 
maximum strain criterion is a strain-based noninteracting 
criterion and, as such, does not differentiate or account 
for possible interactions which may occur between failure 
modes. The Tsai-Wu criterion, on the other hand, can ac- 
count for interactions among the possible failure modes, 
such as an interaction between fiber failure and shear fail- 
ure. The maximum strain criterion is used as the primary 
criterion to determine the failure mode, and the Tsai-Wu 
criterion is used to verify the accuracy of the maximum 
strain criterion and to provide additional insight on possi- 
ble interactions among failure modes. The legend at the top 
of each table defines the notation associated with each 
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Maximum Strain Criterion 


■mil Tsal-Wu Criterion 

Stacking Design Design 1 Failure Load: 27600 Ib/in 

Sequence Number 


(± 45 / 0 / 90)23 1 



1.00 

0.99 





0.83 

0.67 



1.26 

1.10 


( 0/C,) 3 



2.29 
1 . 94 


(±45/Cg)3 


5 



1 . 84 
1.91 


( 02/Cs)3 



2 . 09 
1.93 


(± 45 / 05)3 



1.43 
1 .11 


(±45/03)23 



1.29 
1 . 32 


(±45/C2)23 



1 . 47 
1 . 52 


0.5 1.0 


1.5 2.0 2.5 


Failure Load ( Normalized By Design 1 ) 

( Maximum Strain Criterion ) 


Table 5.1 Normalized Tensile Failure Load 
Case: L = 20 in. , W = 10 in. , D/W = 1/6 


97 





Maximum Strain Criterion jp | | | j| H ||| Tsai-Wu Criterion 

Stacking Design Design 1 Failure Load: 26400 Ib/in 

Sequence Number 


(±45/0/90) 1 



1.00 

0.99 


( Oe). 



0.76 

0.62 


( C.). 



1.18 

1.04 


( 0/C,), 


{±45/Ce), 


4 



2.21 
1 . 87 



1 . 79 
1 . 85 


( Oj/Cg), 



2.03 
1 . 87 


(±45/0,), 



1 . 38 
1.07 


(±45/03)2, 



1.26 

1.29 


(±45/C2)2, 



1 . 44 
1 . 49 


0.5 1.0 


1.5 2.0 2.5 


Failure Load ( Normalized By Design 1 ) 

( Maximum Strain Criterion ) 


Table 5.2 Normalized Tensile Failure Load 
Case: L = 10 in. , W = 10 in. , D/W = 1/6 
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mu Maximum Strain Criterion iiiiiii Tsai-Wu Criterion 

Stacking Design Design 1 Failure Load: 24100 Ib/in 

Sequence Number 


(±45/0/90)2^ 1 



1.00 

0.99 


( 0e)s 



0.68 

0.59 







1.28 

1.11 


( 0/Ct)» 



2.04 
1 . 92 


(±45/C6)„ 



1.71 
1 . 77 


C O 2 /C,), 



1.90 

1.76 


(±45/0g)„ 



1.38 

1.06 


(±45/02)2, 



1.28 

1.31 


(±45/C2) 


2 8 


9 



1.42 
1 . 47 


0.5 1.0 1.5 2.0 2.5 

Failure Load ( Normalized By Design 1 ) 

( Maximum Strain Criterion ) 


Table 5.3 Normalized Tensile Failure Load 
Case: L = 20 in., W = 10 in., D/W = 1/3 
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table. In each table, the failure loads of each design 
(stacking sequence) are normalized by the failure load of 
the straightline quasi-isotropic design (Design 1) for that 
particular geometry based on the maximum strain criterion. 
Quasi-isotropic laminates represent a conventional , conser- 
vative, current day design philosophy. For each geometry, 
if the normalized failure load for a particular design is 
less than 1.00, that design is not as capable of supporting 
a load as the quasi-isotropic design. Conversely, if the 
normalized failure load is greater than 1.00, that particu- 
lar design is predicted to support a higher load at failure 
than the quasi-isotropic design. 

Table 5.1 shows the results for plates with parameters 
L/W = 2 and D/W = 1/6. These plates are rectangular and the 
actual plate dimensions are L = 20 in . , W = 10 in . , and D = 
1.667 in. (length, width, and hole diameter). The maximum 
strain criterion predicts that the Design 1 laminate will 
experience fiber failure in the 0° layers at the net-section 
hole edge at a load of 27600 Ib/in. The Tsai-Wu criterion 
predicts the same result at virtually the same failure load. 

Considering Design 2 next, the unidirectional fiber 
format laminate (0a)s is one that comes to mind when think- 
ing about using fibers most effectively, i.e., aligned with 
the load direction. However. Design 2 experiences shear 
failure in the matrix, at the hole, near the net-section, at 
a load 17% below that of quasi-isotropic Design 1. This 
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failure results in a crack which runs parallel to the fibers 
from the net-section hole edge to the ends of the plate. 
Shear occurs because of the rapid change in the direction of 
stress as the load is transmitted around the hole. In order 
to accomodate the hole, the uniform load on the ends of the 
plate necessarily must become nonuniform as the hole is ap- 
proached. This generates gradients in one stress component 
that in turn produces gradients in the other stress compon- 
ents. Since there are no shear stresses at the ends of the 
plate, shear stresses are generated when approaching the 
hole. Since a unidirectional fiber format is very weak in 
shear, this laminate is design-limited in shear. The Tsai- 
Wu criterion predicts the same result at a failure load 33% 
lower than that of Design 1. 

The choice of a curvilinear fiber format and the stack- 
ing sequence (Cb)s for Design 3, where the fibers are align- 
ed with the principal stress directions at every point in 
the plate, is a logical progression to overcome the inherent 
problems of Design 2 which led to shear failure. Herein "C" 
is used to denote curvilinear layers. As a result, no shear 
stress can occur in the principal material directions. The 
maximum strain criterion predicts that Design 3 is 26% 
stronger than Design 1 (quasi-isotropic) and 43% stronger 
than Design 2 (unidirectional). Failure of the laminate of 
Design 3 is traced to the relative weakness of thQ matrix in 
the direction transverse to the fibers. Design 3 fails due 
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to strain elongation in the fiber direction generating con- 
traction strains in the transverse direction at the net-sec- 
tion near the hole edge (Poisson ratio effect). The materi- 
al simply separates perpendicular to the fiber direction. 
The Tsai-Wu failure criterion predicts the same result at a 
failure load 10% greater than that of Design 1. 

Again, a logical progression in choice of stacking se- 
quence, to overcome the inherent problems encountered with 
Design 3, is that of Design 4, a (0/C7)s laminate, and 
Design 6, a (02/C8)s laminate, where the "0" denotes 
orthogonal. In Design 4, the fibers in the two "0" layers 
are aligned orthogonally everywhere to the 14 load-bearing 
curvilinear layers to prevent failure of the material 
perpendicular to the fiber direction. Similarly, in Design 
6 the fibers in the four "0" layers are aligned orthogonally 
everywhere to the 12 load-bearing curvilinear layers to 
alleviate failure of the material perpendicular to the fi- 
bers at the net-section away from the hole edge. The failure 
mechanism for Design 4 and Design 6 is fiber failure in the 
curvilinear layers at the net-section near the hole edge. 
Design 4 supports the highest fiber direction tensile stress 
in the curvilinear layers among the designs considered, and 
shows that two orthogonal layers are all that is necessary 
to eliminate failure of the material perpendicular to the 
fibers. The maximum strain criterion predicts that Design 4 
is 129% stronger and Design 6 is 109% stronger than Design 
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1. This is a dramatic increase in strength. The Tsai-Wu 
criterion also predicts about a factor of 2 strength in- 
crease relative to Design 1. Unfortunately, any attempt to 
fabricate an orthogonal grid of fibers would probably result 
in less than perfect orthogonal alignment relative to the 
curvilinear fibers. If the alignment was less than perfect 
shear stresses would be 'introduced which would lead to mat- 
rix failure. Even if the orthogonal alignment was perfect, 
if the load was other than pure axial loading (i.e., a shear 
component was present) , shear stresses would be generated 
within the laminate and this would lead to laminate failure 
due to shear. Otherwise, Design 4 and Design 6 are excell- 
ent designs for the case of pure axial load. 

The choice of a more practical design, especially for 
easier fabrication and to anticipate inadvertent introduc- 
tion of shear, is that of Design 5, a (±45/Ce)8 laminate, 
which includes 12 load-bearing curvilinear fiber layers. 
Design 7, a (±45/0e)s laminate, is a comparable laminate 
which includes 12 load-bearing straightline fiber layers 
instead of 12 curvilinear layers. Both designs include 4 
straightline fiber layers at ±450 (two positive and two neg- 
ative to maintain symmetry) which provide the desired resis- 
tance to shear. These layers also help carry axial load. 
Use of the ±45° layers effectively eliminates the problem of 
shear failure in the matrix. The iteration scheme, discuss- 
ed in Chapter III, is used to determine the necessary fiber 
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orientation of the curvilinear layers to maintain fiber 
alignment in the principal stress direction. It is impor- 
tant to note that the directions of the curvilinear fibers 
at a given point in the (±45/Ce)3 laminate are not the same 
as the directions of the curvilinear fibers in the (Ce)s 
laminate, even though the notation "C" is used in both 
cases . ^ 

The maximum strain criterion predicts that Design 5 is 
84% stronger than Design 1, that Design 7 is 43% stronger 
than Design 1, and that Design 5 with its curvilinear fiber 
format has a distinct margin of strength over Design 7 with 
its straightline fiber format. Failure for Design 5 occurs 
in the curvilinear layers, and for Design 7 occurs in the QO 
layers, with fiber failure at the net-section near the hole 
edge as the mechanism for both designs. The Tsai-Wu criter- 
ion predicts the same failure mode for both designs, and 
predicts a 91% strength increase for Design 5. and an 11% 
strength increase for Design 7. both relative to Design 1. 
Moreover, for Design 7, the Tsai-Wu criterion indicates 
interaction between shear failure and fiber failure, with 
fiber failure controlling the response at the hole near the 
net-section. This interaction may account for the increased 
difference in strength (91% versus 11%) between Design 5 
(curvilinear) and Design 7 (straightline). Also, a possible 
difference is that the transmission of the load around the 
hole through the net-section is accomplished differently 
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with the curvilinear fiber format than with the straightline 
fiber format. This will be discussed in more detail later 
in this chapter. 

Care should be exercised not to make comparisons among 
the various designs based on factors which may be too narrow 
and restrictive. All of the designs were kept constant with 
respect to weight by maintaining each design as 16 layer 
laminates. Strength as the only factor to be considered in 
making a best design choice would be too narrow. Based on 
the maximum strain criterion, if Design 6 was to be chosen 
over Design 5 because of its 25% strength margin, the dis- 
advantages associated with possible shear load introduction 
would have been ignored. As can be seen from Table 5.1, 
comparisons based on the Tsai-Wu criterion would lead to the 
same conclusions. However, according to the Tsai-Wu criter- 
ion, possible interactions between fiber failure and shear 
failure favors Design 5. 

Finally, considering Design 8, a (±45/02 )2s laminate, 
and Design 9, a (±45/C2)2a laminate, both criteria predict 
that the curvilinear fiber format will have a greater load 
capacity than the straightline fiber format. Design 8 and 
Design 9 have more ±45° layers, at the expense of less 0° or 
less curvilinear load-bearing layers than Designs 4-7. For 
Design 8. the failure mode is predicted to be fiber failure 
at the net-section near the hole edge in the OO layers. 
According to the maximum strain criterion. Design 8, with 
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less QO load-bearing layers, has a lower failure load margin 
than Design 7, 29% versus 43%, relative to Design 1, the 
baseline quasi-isotropic laminate. However, the Tsai-Wu 
criterion indicates interaction between fiber failure and 
shear failure and predicts that Design 8 will have a higher 
failure load margin than Design 7 (32% versus 11%). Any in- 
dication of possible interaction should not be dismissed 
lightly . 

For Design 9, the failure mode is predicted to be fiber 
failure at the net-section, near the hole edge, in the curv- 
ilinear layers . The maximum strain criterion predicts a 
strength increase of 47% for Design 9 relative to Design 1, 
and a failure load margin of 18% over Design 8, a straight- 
line design. The Tsai-Wu criterion predicts that Design 9 
is 52% stronger than Design 1 and 20% stronger than Design 
8. This agrees with the comparison and trend established 
for Design 5 and Design 7, where the curvilinear fiber 
format has superior load capacity over the similar design 
straightline fiber format. As is the case for Design 5 and 
Design 7, with Design 8 and Design 9, at the net-section 
near the hole edge, the failure mechanism (fiber failure), 
layer orientation, and fiber alignment are virtually iden- 
tical. However, as mentioned above, the curvilinear fiber 
format seems to be able to transmit a higher load around the 
hole and through the net-section than the straightline fiber 
format can . 
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Table 5.2 shows the results for plates with parameters 
L/W = 1 and D/W = 1/6, namely, square plates. The actual 

plate dimensions are L = 10 in., W = 10 in., and D = 1.667 
in. (length, width, and hole diameter). For this geometry, 
the quasi-isotropic design (Design 1) failure load determin- 
ed by the maximum strain criterion is used as a baseline 
value with which to normalize and compare the failure loads 
of the other designs, including the failure loads predicted 
by the Tsai-Wu failure criterion. As shown in Table 5.1 for 
the rectangular plate, the same favorable results for the 
curvilinear designs relative to the straightline designs 
were obtained for the square plate geometry. For example, 
as seen in Table 5.2, the Design 5 square laminate, a 
(±45/C6)a, has a normalized failure load of 1.79 relative to 
the Design 1 laminate, compared to 1.84 for the rectangular 
laminate, and fiber failure is predicted to occur in the 
curvilinear layers at the net-section near the hole edge, 
the problem of shear failure in the matrix being eliminated. 
The square straightline laminate Design 7. a (±45/0e)s, has 
a normalized failure load of 1.30, compared to 1.43 for the 
rectangular laminate, and fiber failure is predicted to 
occur in the Qo layers at the net-section near the hole 
edge. Interaction between shear failure and fiber failure 
causes fiber failure in this laminate at a load 41% below 
that of the curvilinear laminate Design 5. 

Tables 5.3 and 5.4, which consider rectangular and 
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square plates with holes twice as large as the other two 
cases, continue the pattern established in Tables 5.1 and 
5.2. The results show clearly that the curvilinear designs 
are stronger, having higher tensile load capacities relative 
to the straightline designs. Like the other plate geome- 
tries, each curvilinear design exhibits the same relative 
improvement over its s'traightline counterpart. Table 5.3 
shows the results for plates with parameters L/W = 2 and D/W 
= 1/3. The actual plate dimensions are L = 20 in. , W = 10 
in., and D = 3.333 in. (length, width, and hole diameter). 
Failure loads are normalized by that of the quasi-isotropic 
Design 1 laminate. Table 5.4 shows the results for plates 
with parameters L/W = 1 and D/W = 1/3. The actual plate 
dimensions are L = 10 in., W = 10 in., and D = 3.333 in. 
(length, width, and hole diameter). Failure loads are nor- 
malized by that of the quasi-isotropic Design 1 laminate. 
An increase in failure load for the curvilinear designs over 
the straightline designs is realized for all plate geometric 
configurations . 

It is important to note that even though the strengths 
of the curvilinear designs relative to the straightline de- 
signs were basically the same for all four plate geometries 
considered, the absolute strengths of the various geome- 
tries, both straightline and curvilinear, are a function of 
plate geometry. Table 5.5 shows the failure load of all de- 
signs normalized by a single failure load, namely the fail- 
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Table 5.5 


Normalized Tensile Failure Load (Geometry Comparison)^ 


Stacking 

Sequence 

Design 

Number 

L/W = 2 
D/W = 1/6 

L/W = 1 
D/W =1/6 

L/W = 2 
D/W =1/3 

L/W = 1 
D/W = 1/3 

(±45/0/90) 

2a ^ 

1.00 

(0.99) 

0.96 

(0.95) 

0.88 

(0.87) 

0.74 

(0.73) 


2 

0.83 

(0.67) 

0.72 

(0.59) 

0.59 

(0.52) 

0 . 43 
(0.38) 


3 

1.26 

( 1 . 10 ) 

1.13 

( 1 . 00 ) 

1 . 13 
(0.98) 

0.75 

(0.69) 

( O/C 7 ), 

4 

2.29 

(1.94) 

2.11 

(1.79) 

1.79 

( 1 . 68 ) 

1 . 40 
(1.32) 

(±45/Cg)3 

5 

1.84 

(1.91) 

1.71 

(1.77) 

1 . 50 
(1.55) 

1 . 18 
(1.23) 

( Oz/Cg), 

6 

2.09 

(1.93) 

1.94 

(1.79) 

1.67 

(1.54) 

1 . 32 
(1.23) 

(±45/0g). 

7 

1 . 43 
( 1 . 11 ) 

1.32 

(1.03) 

1.21 

(0.93) 

0 . 94 
(0.70) 

(± 45 / 02)23 

8 

1.29 

(1.32) 

1.21 

(1.24) 

1 . 12 
(1.15) 

0.89 

(0.90) 

(± 45 / 02 ) 2 , 

9 

1 . 47 
(1.52) 

1 . 38 
(1.43) 

1 . 25 
(1.29) 

0 . 99 
(1.02) 


I Normalized by Design 1 failure load (maximum strain cri- 
terion), L/W = 2, D/W = 1/6. Failure loads determined by 
the Tsai-Wu failure criterion are specified in parentheses. 
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ure load of the quasi-isotropic rectangular plate with a 
small hole (L/W = 2, D/W = 1/6, Design 1. Table 5.1) as 
predicted by the maximum strain criterion. As can be seen, 
the failure loads going from Tables 5.1 to 5.4 decrease for 
any given design. For example, if the normalized values of 
failure loads for Design 5 (curvilinear fiber format) and 
Design 7 (straightline fiber format) are compared progres- 
sively from Table 5.1 through 5.4, this statement is seen to 
be true. All of the data shown for each geometric configur- 
ation in Tables 5.1 through 5.4 is combined in Table 5.5 to 
illustrate the influence of geometry on failure load. The 
normalized values of the failure loads predicted by applica- 
tion of the Tsai-Wu criterion are enclosed in parentheses. 
Designs which have the largest tensile load capacities are 
realized by laminates with the geometric parameters L/W - 2 
and D/W = 1/6. The geometries in order of decreasing ten- 
sile load capacity are; L/W = 1 and D/W = 1/6; L/W = 2 and 
D/W 3 1/3; L/W = 1 and D/W = 1/3. 

Iteration Procedure And Convergence Of Fiber Directions 

The iteration procedure used to determine the direction 
of the fibers in the curvilinear layers for fiber alignment 
with the principal stress directions was described in some 
detail in Chapter III. In Table 5.6, the fiber directions 
in four representative elements for a laminate with a stack- 
ing sequence of (±45/C6)s and geometric parameters L/W = 1 
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Table 5.6 


Iteration 

Convergence Of Fiber 
Element 28 Element 81 

Directions i 
Element 154 

Element 177 

1 

- 3.74 

-12.60 

0.15 

- 

0.49 

2 

-12.36 

-41.15 

1 . 17 

- 

1 . 69 

3 

-15.65 

-42.30 

1.95 

- 

2 . 34 

4 

-16.91 

-42.93 

2.48 

- 

2.70 

5 

-17.34 

-42.96 

2.82 

- 

2.89 

6 

-17.44 

-42.96 

3.04 

- 

2.99 

7 

-17.41 

-42.96 

3.19 

- 

3 . 05 

8 

-17.35 

-42.99 

3.29 

- 

3 . 08 

9 

-17.29 

-43.01 

3.36 

- 

3.09 

10 

-17.23 

-43 . 04 

3.40 

- 

3. 10 

11 

-17.19 

-43.04 

3.44 

- 

3 . 10 

12 

-17.17 

-43 . 06 

3.46 

- 

3.11 

^ Fiber directions are 

in units of 

degrees 




and D/W = 1/3 are listed as a function of iteration number. 
Location of the four elements in the quarter plate is shown 
in Figure 5.1. Equation 3.7, page 75, is used as the con- 
dition to be satisfied for convergence of the fiber align- 
ment angle. An application of equation 3.7 provides results 
as shown in Table 5.7 where data from Table 5.6 and itera- 
tions 11 and 12 are used to check for convergence. As can 
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Figure 5.1 Location Of Four Elements 
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be seen from the table, convergence is quite rapid, and can 
actually be achieved within 6 or 7 iterations. 


Table 5.7 

Relative Error In Current Principal Stress Directioni 
Element 28 Element 81 Element 154 Element 177 


0.0012 0.0005 0.0058 0.0032 


I Tabular values obtained using data of previous iteration 

II and current iteration 12 from Table 5.6 and applying 
equation 3.7, p. 75. Convergence is assumed as achieved 
when the relation given by equation 3.7 is less than 0.01. 


Stress Continuity Across Element Boundaries 

The question of the accuracy of the method of analysis 
could be a matter of concern. Finite element methods are 
well established methods for providing numerical results. 
The effectiveness of the eight-node isoparametric element in 
general applications to satisfy convergence requirements, 
compatability , and stress continuity is well established in 
the literature. Related discussions will not be included. 
However, in this problem, because of the idealization, the 
fiber orientation varies from element to element, causing 
possible discontinuities in the stress field. It is there- 
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fore logical to ask about the convergence of the results. 
The standard approach is to compute results using a finer 
mesh and compare the results with those of the established 
mesh. If the results from the finer mesh compare well with 
the results of the established mesh, the answers are assumed 
to have converged. Here a finer mesh results in less of a 
jump in fiber orientation from element to element, and 
thereby, perhaps less problems for convergence of the stress 
field. To investigate this point, the standard mesh size 
density was doubled, from 192 elements to 384 elements in a 
quarter plate analysis. The circumferential force resultant 
Na around the hole edge was used as a measure of stress 
field continuity. The 192 element mesh resulted in 12 
elements around the quarter hole, while the 384 element mesh 
resulted in 24 elements around the quarter hole. Data for 
the laminate design case (±45/C6)a, L/W = 1, D/W = 1/3, were 
used for the comparison. In Figure 5.2 the force resultants 
for the two meshes are plotted against element position 
around the hole given by the angle relative to the x-axis . 
The results for the two meshes are virtually superimposed 
and show neglible difference in results. The jump in the 

force resultant from element to element is evidenced by the 
fact that the results are smooth and continuous and is 
negligible. Data used to plot the curves of Figure 5.2 are 
shown in Table 5.8. 
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Figure 5.2 Influence Of Mesh Sine On Results 
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Table 5.8 


Circumferential Force Resultant Data 
Laminate Design Case: (±45/Ce)s L/W = 1 D/W = 1/3 
192 Element Mesh^ 384 Element Meshi 


Element 

Number 

Ns 

(Ib/in) 

Element 

Number 

Ns 

(Ib/in) 

Angle: i 

( Degrees 



1 

- 514.96 

1 . 875 

1 

- 512.32 



3.75 



17 

- 507.70 

5.625 



33 

- 492.50 

9.375 

17 

- 482.16 



11.25 



49 

- 467.94 

13.125 



65 

- 432.21 

16.875 

33 

- 410.63 



18.75 



81 

- 382.92 

20 . 625 



97 

- 317.33 

24.375 

49 

- 280.20 



26 . 25 



113 

- 233.10 

28 . 125 



129 

- 123.39 

31.875 

65 

- 54.29 



33.75 



145 

50.16 

35 . 625 



161 

278.15 

39.375 

81 

417.13 



41 . 25 



177 

546.82 

43 . 125 



193 

829.37 

46 . 875 

97 

975.91 



48.75 



209 

1126.82 

50.625 



225 

1463.28 

54.375 

113 

1643 . 00 



56 . 25 



241 

1813.34 

58 . 125 



257 

2164.73 

61.875 

129 

2336.73 



63 . 75 



273 

2502.58 

65.625 



289 

2811 . 98 

69.375 

145 

2946.92 



71 . 25 



305 

3081 . 38 

73.125 



321 

3301.78 

76.875 

161 

3383.61 



78.75 



337 

3468.43 

80 . 625 



353 

3579.12 

84.375 

177 

3605.42 



86 . 25 



369 

3634.28 

88.125 


I Data for each element taken at Gauss point 5 . 
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A further comparison was made by simply examining the 
circumferential force resultant for the three sets of ele- 
ments bordering the edge of the hole (12 elements on the 
boundary, 12 elements next to those boundary elements, and a 
third set of 12 elements next to the second set) . Essen- 
tially these elements form three concentric rings of ele- 
ments around a quarter of the hole. Data for the laminate 
design case (±45/C6)8, L/W = 1, D/W = 1/3, were also used 
for the comparison. In Figure 5.3 these results for the 
three sets of elements are plotted against element position 
around the hole given as the angle relative to the x-axis . 
The jump in the circumferential force from element to ele- 
ment, both circumferentially and radially as shown by the 
curves, is minimal, as evidenced by the smoothness of each 
curve, and the order among the three curves. Data used to 
plot the curves of Figure 5.3 is shown in Table 5.9. It is 
clear from the results of using two mesh densities, and from 
the smoothness both radially and circumferentially, that the 
original 192 element mesh for the quarter plate is suffi- 
cient. 
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N,: Circumferential Force Resultant (x 100 Ib/in) 

Figure 5.3 Change In Force Element To Element 
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Table 5.9 


Circumferential Force Resultant Datai 
Laminate Design Case: (±45/C6)s L/W = 1 D/W = 1/3 


Element 

Number 

Ns 

(Ib/in) 

Element 

Number 

Ns 

( Ib/in) 

Element 

Number 

Ns 

( Ib/in) 

1 

- 512.32 

2 

- 483.86 

3 

- 454.18 

17 

- 482.16 

18 

- 449.19 

19 

- 415.53 

33 

- 410.63 

34 

- 368.89 

35 

- 328.00 

49 

- 280.20 

50 

- 226.64 

51 

- 172.00 

65 

- 54.29 

66 

15.93 

67 

123.96 

81 

417.13 

82 

463.80 

83 

546 . 98 

97 

975.91 

98 

1025.66 

99 

980.21 

113 

1643.00 

114 

1569.01 

115 

1475 . 45 

129 

2336.73 

130 

2174.00 

131 

2023.99 

145 

2946.92 

146 

2735.79 

147 

2544.69 

161 

3383.61 

162 

3164.88 

163 

2956 . 41 

177 

3605 . 42 

178 

3393 . 46 

179 

3183.57 

1 Data for each element taken 

at Gauss 

point 5 . 
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Influence Of The Curvilinear Fiber Format On The Load 
Carrying Capacity 

The curvilinear design (±45/Cs)s and the straightline 
design (±45/Oe)s look very similar near the net-section hole 
edge. The fibers in the curvilinear layers of the curvilin- 
ear design pass by the net-section perpendicular to a line 
from the hole edge to the plate edge. For the straightline 
design, the fibers in the 0® layers are also perpendicular 
to a line from the hole edge to the plate edge. Both de- 
signs have ±450 layers with identical orientation at the 
net-section hole edge. Locally, near the net-section hole 
edge both designs are predicted to have fiber failure at 
this point. However, the curvilinear design is able to 
carry a much higher tensile load than the straightline 
design. Why is this so? 

Figures 5.4 and 5.5 are presented in order to provide 
insight as to why the curvilinear design has a larger ten- 
sile load capacity than the straightline design. Contour 
lines of stress in the x-direction are shown in Figure 5.4 
for the curvilinear layers in the design (±45/C6)s. Contour 
lines of stress in the x-direction are shown in Figure 5.5 
for the 00 layers in the design (±45/0e)a. The curvilinear 
design has fewer contour lines grouped near the net-section 
hole edge, and the contour lines are generally spread far- 
ther apart. This behavior indicates a more gradual change 
in the stress gradients at the net-section. This can be 
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interpreted to mean that the stress concentration is lower 
near the net-section hole edge for the curvilinear design 
than it is for the straightline design. Figure 5.4 shows 
that the stress near the net-section hole edge is approxi- 
mately 3750 psi due to a tensile load of 7200 lb for the 
curvilinear design. Figure 5.5 shows that the stress near 
the net-section hole edge is approximately 4250 psi due to a 
load of 7200 lb for the straightline design. The plot of 
contour lines in Figure 5.4 indicates that the curvilinear 
design carries the load along stress trajectories in a 
streamline pattern which flow around the hole. 

An even better example of this is shown in the follow- 
ing two figures. Figures 5.6 and 5.7 show plots of contour 
lines of stress in the x-direction in the curvilinear layers 
for the orthogonal curvilinear designs (02/C6)s and (0/C7)s. 
Here it is clear the orthogonal layers help to move the 
stress away from the hole edge, reducing the gradients 
around the hole and allowing the load to be transmitted more 
smoothly around the hole. Though it is not needed, because 
the (0/C7)a laminate is sufficient, the value of doubling 
the orthogonal layers is clear. Those layers help move the 
load away from the hole. 

This concludes the chapter on tensile loading. It is 
clear that the variable fiber orientation is beneficial. 
The next chapter addresses the influence of these designs on 
the buckling loads. 
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Fiffure 5.6 Plot. Of Contour Lines Of Stress c% 
For Curvilinear Layers 
Case: (02/Cs)3 . L/W = 1 D/W = 1/3 
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Figure 5.7 Plot Of Contour Lines Of Stress ax 


For Curvilinear Layers 
Case: (0/C7)s L/W = l D/W = 1/3 
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CHAPTER VI 


BUCKLING RESISTANCE OF THE PLATE DESIGNS 

Establishment of criteria applicable to the design and 
analysis of plates subjected to tensile loads was of primary 
concern in the previous chapter. Having considered and com- 
pared plate designs which carry tensile loads, it is of in- 
terest to determine the buckling resistance of these de- 
signs. While it may be possible to design the curvilinear 
fiber format to best resist buckling, that is not the intent 
here. Here, simply, the tensile load designs were checked 
for their resistance to buckling. The buckling analysis was 
conducted to determine the extent to which the tensile-load- 
designed curvilinear fiber format could influence and per- 
haps improve the capacity of the plate to resist buckling. 

The plates were subjected to in-plane compression loads 
and the critical compression load was determined. The crit- 
ical compression load of each plate was defined to be that 
load which would cause the plate to buckle into its first 
buckling mode shape. Figure 6.1 shows a typical first buck- 
ling mode for square plates, and Figure 6.2 shows a typical 
first buckling mode for rectangular plates with length to 
width ratios of 2.0. For the square plate the buckling mode 
is a single half-wave in both directions, while for the rec- 
tangular plate there are two half-waves in the lengthwise 
direction . 
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ORIGINAL PAGE IS 
OF POOR QUALITY 


Figure 6.1 


First Buckling Mode Shape 
Square Plate 
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Figure 6.2 First Buckling Mode Shape 
Rectangular Plate: L/W = 2.0 


12S 


Tables 6.1 through 6.4 present the critical buckling 
loads of 36 different plate designs. There are four differ- 
ent geometries, each having the nine stacking sequences con- 
sidered in Chapter V. The critical buckling loads of the 
nine stacking sequences for each geometry are normalized by 
the critical buckling load of the quasi-isotropic design of 
each geometry. The buckling analysis results for the square 
plates were verified by comparing the results of the (Os) a 
laminates in Tables 6.2 and 6.4 with the results obtained by 
Nemeth [8]. The buckling analysis results for the rectangu- 
lar plates in Tables 6.1 and 6.3 were verified by running a 
test case on a rectangular plate with a hole diameter to 
plate width ratio equal to 0 . 4 with a stacking sequence of 
(±458 )s and then comparing those results with the results 
obtained by Nemeth [8]. In both cases the results obtained 
in this study agreed very well with the results given in 
Reference 8 . 

The results given in Tables 6.1 through 6.4 show that 
the curvilinear fiber format is as good as a straightline 
format as regards buckling resistance. Interestingly 
enough, none of the designs considered, straightline or 
curvilinear, is better than the quasi-isotropic design in 
resisting buckling. For example, considering the data of 
Table 6.1, the (08)s design resists only 38% of the buckling 
load of the quasi-isotropic design. Allowing the straight 
fibers to become curvilinear does little to improve the 
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Stacking Desi 
Sequence Nuin‘ 

(±45/0/90)2, 1 
( Oa)a 2 
( Ce), 3 
( 0 /C,), 4 
(±45/Cg), 5 

( Oa/Cfi), 6 
(±45/0^), 7 
(±45/02)2, 8 
(±45/C2)2, 9 

Ta 


Case : 














Design 1 Buckling Load: 1443.33 lb. 


Stacking Design 
Sequence Number 


(±45/0/90)2, 1 






( 0/C,), 4 


(±45/Cg), 5 


( O2/C6), 6 


(±45/0g), 7 


(±45/02)2, 8 


(±45/C,) 


2 ^ 2 a 


0.0 


1 . 00 


0.63 


0.66 


0.66 


0.94 


0.64 


0.93 


0.98 


0.99 


0.5 


1.0 


1.5 


Buckling Load 
( Normalized By Design 1 ) 


2.0 


Table 6.2 Critical Buckling Load 
Case: L = 10 in., W = 10 in., D/W = 1/6 
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Design 1 Buckling Load: 1727.94 lb. 

Stacking Design 
Sequence Number 


(±45/0/90)2. 1 


t Oa). 


( Ca). 


( 0/C,). 4 


(±45/Cg). 5 


( O 2 /C,). 


(±45/0e). 7 


(±45/03)2. 8 


(±45/C2)2„ 


0.0 


1.00 


0.24 


0.32 


0.72 


0.93 


0.62 


0.87 


0.94 


0.94 


0.5 


1 . 0 


1 . 5 


Buckling Load 
( Normalized By Design 1 ) 


2 . 0 


Table 6.3 Critical Buckling Load 
Case: L = 20 in., W = 10 in.. D/W = 1/3 
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load. Adding orthogonal layers to make a (0/C7)s laminate 
greatly improves the buckling capacity relative to the (Oa)s 
and the (C8)s designs, but it is not as good as the quasi- 
isotropic laminate. Adding the ±45® layers to form the 
(±45/Cs)s laminate dramatically improves the buckling resis- 
tance. with the ±450 layers having a substantial influence. 
The (02/C6)a laminate is very similar to the (0/C7)s lami- 
nate. As a comparison with the (±45/Cb)s laminate, the 
(±45/06 )s laminate shows very similar buckling capacity. 

> the (±45/02 )28 laminate and the (±45/C2)2s laminate 
show similar characteristics and are both as good as the 
quasi-isotropic designs. Aside from understanding the 
influence of the curvilinear layers, the results in Table 
6.1, and similarly, the other three tables, indicate that 
the ±450 layers strongly influence buckling resistance. 

By examining the four tables it is also clear that 
there is a geometric influence on the normalized results. 
For example, for the rectangular plate with the small hole 
(Table 6.1), the (C8)s laminate has only 39% of the buckling 
capacity of its quasi-isotropic counterpart. For a square 
plate with the small hole (Table 6.2), the (C8)s laminate 
has 66% the buckling capacity. While this laminate may not 
be of practical interest, it does demonstrate the influence 
of geometry. 

A,t this point the tensile and buckling capacities of 
particular plates can be compared to provide insight into 
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the overall gains of using the curvilinear fiber format. 
The tensile and buckling capacities of several laminates and 
the range of geometries are compared in Table 6 . 5 and Table 
6.6. As before, the buckling capacity of the quasi- 
isotropic laminate in each case is defined to be unity. 
From the table it is clear that the curvilinear design has 
an advantage. While the buckling capacity of both the 
(±45/06) a laminate and the (±45/Ce)s laminate are about 5% 
less than the buckling capacity of the quasi-isotropic lami- 
nate as listed in Table 6.5, and even though both laminates 
carry more in tension than the quasi-isotropic plate, the 
(±45 /Cg)b laminate is better in tension than the straight- 
line counterpart as listed and shown in Table 6.6. With the 
buckling loads being the same, the increased capacity of the 
curvilinear design makes it superior. Similar comments can 
be made regarding the (±45/02)28 laminate and the (±45/02) 2 s 
laminate . 

In summary, it can be said that the curvilinear fiber 
format used in the curvilinear designs do not degrade the 
buckling performance. 
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Table 6.5 


Critical Buckling Loadi 



L/W = 2 
D/W = 1/6 

L/W = 1 
D/W = 1/6 

L/W = 2 
D/W = 1/3 

L/W = 1 
D/W = 1/3 

(±45/0/90)28 

1.0 

1.0 

1.0 

1.0 

(±45/06)8 

0.95 

0.93 

0.87 

0.84 

(±45/C6)s 

0.96 

0.94 

0.93 

0.87 

( ±45/02 ) 2 8 

1.00 

0.98 

0.94 

0.90 

(±45/C2 ) 28 

1.00 

0.99 

0.94 

0.93 



Tensile 
L/W = 2 
D/W = 1/6 

Table 6.6 
Load Capacity 
L/W = 1 
D/W = 1/6 

■1 

L/W = 2 
D/W = 1/3 

L/W = 1 
D/W = 1/3 

(±45/0/90)28 

1 . 0 

1 . 0 

1.0 

1 . 0 

(±45/06)3 

1 . 43 

1 . 38 

1.38 

1.27 

(±45/C6)8 

1.84 

1.79 

1.71 

1.60 

(±45/02 ) 2 8 

1.29 

1.26 

1.28 

1.20 

(±45/02)28 

1 . 47 

1 . 44 

1.42 

1.33 

1 Designs are 

normalized 

by the quasi- 

isotropic 

design for 


each geometry 
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CHAPTER VII 


SUMMARY. CONCLUSIONS AND RECOMMENDATIONS 


SuiDiDary 

This study has addressed and investigated the issue of 
using a curvilinear fiber format to increase the structural 
efficiency of layered fiber-reinforced, simply supported, 
flat plates with centrally located circular holes. The is- 
sue of how to determine the orientation and alignment of the 
fibers from point to point in a laminated plate structure , 
and how to evaluate the gains and improvement in using the 
format, have been addressed. 

In-plane tensile loading was applied to two opposing 
plate edges and the fiber directions of the curvilinear lay- 
ers were determined. Specifically, an iteration scheme and 
procedure, in conjunction with a finite element discretiza- 
tion of the plate, was formulated and used to determine the 
fiber orientation and alignment in some or all of the layers 
such that in those layers computed values of the principal 
stress directions and the principal material directions co- 
incided. Evidence of the convergence of the scheme was also 
presented . 

In establishing a basis for comparison, a maximum 
strain failure criterion and the Tsai-Wu failure criterion 
were both applied and used to predict and determine the max- 
imum in-plane tensile loads which could be carried and the 
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failure modes of the various designs. The Tsai-Wu failure 
criterion was applied to support, check and compare the 
numerical results obtained with the maximum strain failure 
criterion. The Tsai-Wu failure criterion agreed fully and 
provided some additional insight on the interactive nature 
of some failure mechanisms. Increase in capacity was pre- 
dicted by both failure criteria. 

Various geometries were examined and results were com- 
pared to a quasi-isotropic laminate with a similar geometry. 
A quasi-isotropic laminate was chosen as a basis since it 
represents a rather conventional and often-used design. 
Several idealistic designs were also considered, e.g., a 
( 0 /C 7 )a laminate, and these were found to show about a fac- 
tor of two increase in tensile load capacity. Due to possi- 
ble manufacturing difficulties and the inability to resist 
shear loading, variants of these ideal designs were also 
considered. This involved replacing the orthogonal layers 
with ±450 layers. Those laminates have straightline 
counterparts which are used in current design philosophy. 
Compared to those straightline counterparts, and the quasi- 
isotropic basis, the curvilinear designs showed improved 
tensile load capacity. 

It was found that neither the hole size or the aspect 
ratio of the plate had any correlative effect on increasing 
the tensile strength of the curvilinear designs over that of 
the straightline designs. Numerical results obtained on in- 
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plane tensile strength of the curvilinear designs and 
straightline designs were highest for the plate parameter 
case D/W = 1/6 and L/W = 2, followed in descending order by 
the cases D/W = 1/6 and L/W = 1, D/W = 1/3 and L/W = 2, and 
D/W = 1/3 and L/W = 1 . Contour plots of the force resultant 
in the x direction indicated that a possible reason for the 
improved capacity of the curvilinear design is that the load 
path is directed away from the hole edge, with the curvilin- 
ear fibers causing the load to 'flow' around the hole. 

Though tensile capacity is important, compression load 
capacity, particularly in the sense of resisting buckling. 
Is equally important. Thus the designs which evolved for 
improving tensile loading were checked with regard to their 
resistance to buckling. Again, buckling resistance was com- 
pared to a quasi-isotropic basis. Numerical results for 
both the curvilinear and straightline designs were obtained 
by applying a finite element method in the commercial code 
Engineering Analysis Language to establish a basis for com- 
paring relative resistance to buckling when the designs un- 
der simply supported edge conditions were subjected to uni- 
axial compressive loading on two opposing edges. This por- 
tion of the study indicated that the curvilinear designs had 
little influence on buckling loads compared to their 
straightline counterparts. Interestingly enough, no design 
was better than the quasi-isotropic design in resisting 
buckling, and the presence of ±45° layers improved the buck- 
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ling resistance for all of the designs. 


Conclusions 

Of primary importance, it is concluded that use of a 
curvilinear fiber format improves tensile load capacity and 
does not degrade buckling capacity. 

Of secondary importance, it is concluded that 

1. The iteration scheme for finding fiber orientation 
worked well and should work for other structural elements, 
o.g., curved panels. 

2. In general, the maximum strain and the Tsai-Wu 
failure criteria predicted identical load capacities in 
tension. When there were differences it was due to 
prediction of interaction effects with the Tsai-Wu 
criterion . 

3. The addition of and/or presence of ±45° layers 
enhances buckling resistance. 

Recommendations - Suggestions For Further Study 

A natural departure for further study is to consider 
designs which include nonsymmetrical stacking sequences. 
The area of interest would be to determine if laminate 
designs which use a curvilinear fiber format could carry 
loads more efficiently than similar designs which use a 
straightline fiber format. Stiffness terras which account 
for coupling effects resulting in twist due to extension 
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force (or extension due to moments) would necessarily 
increase the complexity of the problem. 

Another point of departure for further study would be 
concerned with curved plates (cylindrical sections) and 
shells subjected to in-plane loads as well as transverse 
distributed loads . An area of interest to be investigated 
would involve determining whether and when effects of trans- 
verse shear forces should be considered rather than be 
neglected . 

In this study, design of the curvilinear fiber format 
laminates was based on the application of an in-plane ten- 
sile loading. The fiber directions in the curvilinear layers 
were oriented to provide the most efficient resistance for 
that type of loading. Buckling resistance of the laminates 
was shown to be lower than their levels of tensile strength. 
This establishes response to buckling as a critical design 
mode and buckling strength as a resistance level which the 
laminate must support without failure. Thus, further study 
may be beneficial in generating laminate design criteria by 
applying in-plane compressive loading to determine whether a 
basis could be established to justify that fibers should be 
aligned in a specified direction to increase resistance to 
buckling and have the design remain capable of satisfying 
design requirements for tensile loads. 
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Finally, consideration should be given to fabricating 
specimens to verify the findings of the study. Baseline 
designs, curvilinear designs, and their straightline coun- 
terparts should be fabricated and tested in tension and in 
compression . 
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